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Abstract 


A  comprehensive  research  program  is  under  way  for 
the  characterization  of  electronic  transport  and  dielectr^ 

,  VHnwn  nrcoet r'es  of  thin  insulating  films  cn  Si  which 
™  ^SLScTcIu,  l^orc-nc .  namely  **2.  ^3.  “X*  “* 

their  composites.  The  main  experimental  «*ppro«»che.  a 

i)  Non-destructive  breakdown  of  un-metallized  films 
exposed  rn  B  ccrcna  discharge  in  s  suitable  g-s  (dry 
air,  Nj,  rare,  gas,  ecc.) 

ii)  Locally  destructive,  self-quenching  breakdown  of 
MIS  capacitor  structures. 

lii)  Optical  charge  -  optical  and  thermal  discharge 
studies  of  electronic  traps  in  the  films. 

iv)  Kilovolt  electron-beam  probe  studies  of  the  films, 
including  damage  induced  by  the  beam. 

Accompanying  theoretical  studies  initially  are 
concerned  with  Monte  Carle  calculations  of  hot-electron 
distributions  indued  by  high  electric  fields  in  the  film. 

Very  substantial  progress  was  made  on  this  problem, 

including: 

i)  the  identification  of  useful  staling  laws  for  constant 
mean- free -path  problems,  including  anisotropy  of 
scattering, 

ii)  demonstration  of  the  usefulness  or  one-dimensional, 
rando  ,-walk  simulations  of  three-dimensional  problems, 

and 

ill)  development  of  a  powerful  graphical  technique  for 
studying  distribution  instability  (runaway)  induced  y 
an  energy-dependent  mean  free  path. 
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1  Introduction 


A  comprehensive  research  pregram  is  under  way  for 
tie  characterization  of  electronic  transport  and  dielectric 
freak  down  properties  of  thin  insulating  films  on  Si  which  are 
technologically  important  (SiOg,  Al^,,  Si  and  their 
composites .  Completely  non-destructive  breakdown  of  the  film 
is  readily  achieved  by  charging  the  exposed  l  unmetallized) 
surface  through  either  an  impacting  ion  beam  or  a  corona 
discharge  in  a  suitable  gas.  We  are  concentrating  on  the 
corona  discharge,  which  has  two  strong  advantages:  1)  either 
sign  cl  charging  is  equaxly  convenient,  and  ii)  there  is 
negligible  kinetic  energy  of  impact.  A  direct  comparison 
technique,  employing  a  dummy  metal  sample,  has  been  developed 
hich  makes  it  possible,  for  the  first  time,  to  measure  the 
surface  potential  of  the  exposed  film  during  steady-state 
charging.  A  special  environmental  chamber  with  desiccation 
capability  and  with  gloves  for  manipulation  of  objects  inside 
the  chamber  from  the  outside  has  been  constructed  for  the 
systematic  study  of  corona-induced  breakdown  under  highly 
controlled  conditions  in  a  variety  of  suitable  gases.  A 
second  technique  for  studying  film  breakdown  employs  thin 
metallization,  leading  to  local  removal  of  gate  material, 
on  the  scale  of  microns,  at  the  site  of  a  breakdown,  ending 
m  an  open  circuit  at  that  site.  Thus,  the  breakdown  is  self¬ 
quenching  Initially,  we  are  studying  pj/rolytically  deposited 
A1?03  films,  of  approximately  1000A  thickness,  with  gold 
metallization.  A  third  technique  for  studying  transport 
properties  of  the  film  is  the  charge-discharge  technique.  The 
charging  of  traps  in  the  oxide  film  is  done  optically,  either 
thru  electron-hole  pair  generation  throughout  the  film  bulk 
bv  soft  X-rays  or  by  internal  phot oemiss ion  of  electrons  (or 
holes)  into  the  film  from  either  the  Si  substrate  or  the  metal 
gate.  The  discharging  of  the  traps  is  done  either  by  heat  (at 
elevated  temperatures)  or  by  sub-bandgap  light.  Initially, 
we  are  using  these  charge-discharge  techniques  to  study  thermally 
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grown  SiO„  films,  both  without  further  treatment  and  with  ion 
implants  (Na  and  A1 )  .  The  implants  have  been  carried  out  to 
cent roliab ly  dope  the  oxide  with  electron  traps  and  our  interest 
1E  accordingly  focussed  on  the  trapping  properties  of  the  films. 

A  fourth  technique  for  styding  electronic  transport  through  the 
films  utilises  direct  electron  injection  into  the  film  by  a 
moderate-energy  electron  beam  .energy  in  the  low  kilovolt 
range).  The  beam  probes  the  as-grown  film  and  also  probably 
produces  characteristic  forms  of  damage  (e.g.  at  the  oxide- 
substrate  interface,  and  possibly  throughout  its  range). 

The  latter  will  influence  the  I-V  characteristic  of  the  film  and 
can  also  be  monitored  by  other,  more  standard  techniques. 

Results  have  been  obtained  with  all  four  techniques  and  they  are 
presented  in  the  appropriate  sections.  Since  all  of  the 
experimental  results  are  preliminary  at  this  stage  they  are 
net  highlighted  in  the  present  Introduction. 

The  above  experimental  programs  are  designed  to 
explore  the  film  bulk  and  film-substrate  interface  influences 
on  transport  and  breakdown  These  programs  are  complemented 
by  a  program  of  film  surface  studies  based  on  secondary  electron 
spectrometry.  This  spectrometry  makes  it  possible  to  identify 
atomic  species  on  the  film  surface,  to  explore  the  chemical 
environment  of  a  species,  to  map  out  valence-band  structure 
for  the  host  film  in  the  neighborhood  of  the  exposed  surface, 
and  to  determine  near-surface  plasmon  and  phonon  energies.  A 
specially  designed  spectrometer  has  beer:  built  at  Princeton 
with  both  unusually  high  sensitivity  (better  than  10  monolayer 
-under  optimum  conditions)  and  excellent  spatiai  resolution 
limited  only  by  the  incident  electron  beam  spot  size. 

The  experimental  programs  are  complemented  by  a 
theoretical  program  aimed  particularly  toward  elucidating  the 
physical  mechanisms  'underlying  film  breakdown  and  the  develop¬ 
ment  of  simple  models  for  the  dynamical  spatial  and  temporal 
evolution  of  breakdown  events.  Initially  we  have  made  extensive 
Monte  Carlo  calculations  of  the  behavior  of  hot  electrons  at 


high  fields  in  the  insulating  films.  We  have  shown  that  with  a 
constant  (energy-independent)  mean  free  path  for  energy-loss 
scattering,  there  is  always  a  stable  hot-electron  distribution, 
no  matter  how  anisotropic  the  scattering  is,  short  of  pure 
(100%)  forward  scattering  We  have  developed  scaling  laws  for 
the  average  energy  of  the  stable  distribution  and  for  the 
spatial  development  distance  of  the  distribution.  We  have 
also  successfully  simulated  the  full  three-dimensional  problem 
with  highly  simplified  one-dimens ionai  analogs.  Finally,  we 
have  shown  that  with  an  energy -dependent  mean  free  path  the 
problem  of  stability  (non-runaway)  of  the  distribution  can  be 
successfully  analyzed  by  a  very  simple  graphical  method. 

2.  Non-Destructive  Breakdown  of  Un-Metallized  Films  (Z.  Weinberg 
collaborating) 

A.  Background  Discussion 

The  study  of  charge  transport  in  insulating  films  is 
confronted  with  the  problems  of  achieving  controlled  injection 
of  charge  and  avoiding  destructive  breakdown  of  the  sample. 

One  solution  to  these  problems  which  is  being  vigorously 
pursued  utilizes  atomic  cr  molecular  ions  to  make  direct 
electrical  contact  to  an  exposed,  unmetallized  surface  of  the 
insulator.  Either  a  directed  beam  of  ionized  atoms  may  be 
used,  in  which  case  a  substantial  ion  kinetic  energy  must  be 
reckoned  with,  or  alternately  a  gas  discharge  may  be  used,  in 
which  case  the  ion  kinetic  energy  will  be  negligible,  A 
negative  gas  ion  that  has  an  electron  affinity  smaller  than  that 
of  the  insulator  (which  is  on  the  order  of  1  ev  for  the 
insulators  of  interest)  can,  in  principle,  deposit  its  excess 
electron  directly  into  the  ccnducTicn  band  of  the  insulator. 

If  electron  trapping,  at  the  exposed  surface  or  in  the  film 
bulk, is  not  too  strong,  the  electron  so  injected  can  transport 
through  the  film  to  the  silicon  substrate  where  it  will  be 
collected.  If  the  electron  so  injected  is  trapped  at  or  near 
the  exposed  surface,  this  surface  will  become  charged.  The 
potential  of  the  surface,  during  charging,  can  be  monitored 
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continuously  by  means  of  a  direct  comparison  technique  which 
is  describee  below.  Immediately  following  cessation  of  charging 
the  potential  of  the  charged,  exposed  surface  can  be  monitored 
with  a  vibrating  Kelvin  probe.  Film  breakdown  induced  by 
charging  of  the  exposed  surface  is  completely  non-destructive 
because  of  the  absence  of  lateral  conductivity  along  the  exposed 
surface;  thus,  weak  spots,  which  are  normally  present  in  any 
capacitor  structure,  cannot  discharge  the  electrostatic  energy 
stored  in  the  structure.  Note  that  if  a  negative  ion  with 
electron  affinity  exceeding  that  of  the  insulator  arrives  at 
the  exposed  surface,  the  ion  will  deposit  its  excess  electron 
directly  into  a  localized  surface  state  of  the  insulator,  which 
may  be  the  gas  ion  itself  adhering  to  the  exposed  insulator 
surface  This  situation  is  virtually  equivalent  to  the 
above-discussed  on£ ,  in  which  the  electron  is  first  injected 
into  the  conductioq.  band  of  the  insulator  and  then  captured 
into  a  localized  state  at  or  near  the  exposed  surface.  In 
fact,  it  is  very  difficult  to  distinguish  these  two  situations 
experimentally . 

The  situation  with  regard  to  the  arrival  of  positive 
gas  ions  at  the  exposed  insulator  surface  is  complementary  to 
that  discussed  above  for  negative  ions.  A  positive  gas  ion 
which  has  a  (neutral)  ionization  energy  (as  modified  b$r  iihage 
forces)  that  exceeds  the  sum  of  the  electron  affinity  of  the 
insulator  and  its  band  gap  will  offer  an  empty  state  (below 
the  edge  of  the  valence  band  of  the  insulator)  to  which  .an 
electron  from  the  valence  band  can  make  a  transition;  thus  a 
hole  is  injected  into  the  valence  band  of  the  insulator,  where 
it  may  be  either  trapped  or  mobile.  If  the  hole  is  trapped  at 
or  near  the  exposed  surface,  this  surface  will  become  charged 
and  its  potential  will  rise.  If  the  hole  is  mobile,  it  can 
be  transported  through  the  insulator  into  the  back  contact. 

Even  after  nearly  a  decade  of  MOS  technology,  elementary 
information'*'  beariflg  on  hole  mobility  in  the  valance  band  of 
Si0’2,  including  th$  possibility  of  self-trapping  of  the  hole  in 


non-bonding  oxygen  orbitals,  is  st ill  lacking.  If  the 
(neutral)  ionization  energy  of  ■'■he  incident,  positive  gas  ion 
is  smaller  than  energetically  required  for  hole  injection,  the 
positive  charge  of  the  ion  will  be  deposited  directly  into  a 
localized  surface  stat*-  of  the  insulator,  which  (as  in  the  neg¬ 
ative  ion  case)  may  be  the  gas  ion  itself  adhering  to  the 
exposed  surface  As  before,  the  surface  potential  can  be 
monitored  both  during  and  after  charging. 

Two  types  of  experiments  have  been  initiated-  One 
of  these  employs  a  beam  of  positive  noble  gas  ions  in  a  vacuum. 
The  other  uses  ions  formed  in  a  corona  discharge,  of  either 
positive  or  negative  polarity,  m  a  gas  of  the  chosen  type  at 
approximately  atmospheric  pressure. 

B.  Description  of  Experiments 

For  the  noble-gas  ion-beam  experiments  we  have  designed 
and  constructed  the  apparatus  shown  schematically  in  Fig.  1.  A 
plasma  is  produced  in  a  noble  gas  by  a  low-voltage  arc  dis¬ 
charge  between  a  hot  tungsten  filament  and  an  anode.  The 
arc  is  confined  by  a  short,  narrow  canal  with  metallic  walls, 

and  positive  ions  are  extracted  from  the  plasma  through  an 

3 

aperture  in  the  wall.  The  ions  are  focused  into  a  beam  by 
a  set  of  cylindrical  lenses  The  specimen  chamber  is  evacuated 
to  the  range  of  10  -  10"  t-orr  by  a  turbo-molecular  pump. 

The  noble  gas  is  bled  into  the  plasma  chamber  to  provide  a 
pressure  at  which  the  mean  free  path  is  suitably  small  for 
maintenance  of  the  arc  discharge  A  rotating  sample  holder  is 
capable  of  carrying  several  samples  and  is  equipped  with  a 
feedback-controlled  heater.  We  are  installing  a  vibrating 
Kelvin  probe  for  the  measurement  of  the  surface  potential  of 
the  samples  (as  shown  schematically  at  the  bottom  of  Fig.  1). 

With  this  equipment  we  can  observe,  non destructively ,  the 
current- voltage  characteristics  of  insulating  films  and  the 
dynamics  of  surface  charging  and  decay. 

For  the  crrcna  experiments  we  have  constructed  the 
equipment  shown  schematically  in  Fig  2a.  A  high  voltage,  of 
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either  positive  or  negative  poiarity,  is  applied  to  a  needle 
ej.e:trcde  and  produces  a  corona  discharge  in  a  selected  gas 
at  approximaT ely  atmospheric  pressure  Ions  of  the  corres¬ 
ponding  sign  are  drifted  to  the  sample  by  the  electric  field. 
Inasmuch  as  the  mean  free  path  is  of  the  order  of  10  ^  cm, 
the  ions  are  continuously  thermalized  and  strike  the  sample 
with  only  the  very  small  kinetic  energy  corresponding  to  their 
temperature  The  negligible  energy  of  impact  avoids  a  host 

4 

of  complications  which  are  Just  beginning  to  be  appreciated, 
an  1  constitutes  a  marked  advantage  for  the  corona  deposition 
technique 

The  surface  potential  of  the  sample  under  steady-state 
current-flow  conditions  can  be  determined  by  a  direct  comparison 
scheme  m  the  manner  shown  in  Fig.  2b.  The  sample  and  a  test 
plate  are  arranged  in  symmetrical  manner  When  tl.j  surface 
potentials  of  the  sample  and  test  plate  are  unequal,  the  asymmetry 
introduced  into  the  electric  field  distribution  will  cause  the 
drift  of  ions  to  be  unequal  to  the  two  electrodes  and  an 
unbalance  in  current  will  result-  The  currents  are  balanced 
by  adjusting  the  potential  of  one  of  the  two  electrodes,  whereupon 
the  potential  of  the  sample  and  test  plate  can  be  read  directly 
from  the  test  plate.  This  technique  for  measuring  surface 
potential  during  steady-state  (or  quasi-steady-state )  current 
flow  has  been  successfully  tested  in  our  Device  Physics 
Laboratory  on  such  widely  different  insulators  as  Si02  films 
and  thermoplastic-resin  films  We  feel  that  it  represents 
an  important  advance  in  experimental  technique  in  the  study  of 
the  charging  of  insulating  films,  A  report  on  the  technique 
is  being  prepared  for  publication. 

The  Kelvin-probe  method  of  measuring  surface  potential, 
shown  in  Fig,  2c,  is  considerably  more  sensitive  than  the 
balancing  method  and  is  particularly  useful  for  measuring  the 
discharge  of  the  surface  after  the  ion  source  is  turned  off. 

A  well-known  and  convenient  technique  for  measuring 
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MOS  flatband  voltage  Involves  the  us=  of  a  mercury  probe,  as 
illustrated  in  Fig.  2d.  A  bead  of  mercury  projects  from  a 
glass  capillary  tube  which  is  carried  or  micromanipulators,  and 
this  is  used  to  provide  a  temporary  metallic  ccntact  to  the 
surface  of  the  insulator  without  the  need  for  conventional 
evaporated  metallic  field  plates  which  would  have  to  be 
removed  and  replaced  various  times  during  the  course  of  the 
experiments.  The  flatband  voltage  directly  gives  the  total 
charge  stored  in  the  ir.otu.atcr  in  the  bulk,  and  near  the  silicon 
interface , 

C .  Experimental  Results 

At  the  present  time  our  efforts  are  being  concentrated 

cn  the  corona  discharge  studies  At  the  outset,  several 

interesting  results  have  been  obtained-  We  have  found  that  when 

anSiO^  film,  thermally  grown  cn  a  silicon  substrate,  is  charged 

by  corona  in  air,  the  field  strength  at  breakdown  depends 

on  the  polarity  of  the  discharge,  and  is  approximately 

7  x  10  V/cm  for  positive  corona  and  1-3  x  10^  V/'cm  for  negative 

corona.  These  results  are  m  agreement  with  independen 

measurements  recently  made  at  RCA  Laboratories  ^  Further,  our 

measurements  of  MOS  flatband  voltage  show  that  after  charging 

with  negative  corona,  the  SiO_  has  a  positive  accumulation  of 

^  12  2 

charge  amounting  to  approximately  5  x  10  charges/ cm  ,  such 
as  would  occur  if  holes  were  injected  from  the  silicon  into 
the  valence  band  of  the  insulator  and  were  trapped,  (Rote 
that  this  result  does  not  establish  tbeidentity  of  the  dominant 
current-carrier  in  the  oxide).  Charging  by  positive  corona  in 
air  results  in  a  high-frequency  (1  MHz)  C-V  curve  that  shows 
an  unusually  slew  transition  from  accumulation  to  inversion , 

This  effect  has  been  attributed  in  the  past  to  the  presence 
of  interface  states,^  but  recent  evidence  indicates  that  it 

may  be  the  result  of  lateral  inhomogeneities  in  a  distribution 

.  .  T 

of  impurity  ions  at  the  insulator-semiconductor  interface, 

these  ions  perhaps  having  been  driven  inward  from  the  free 
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surface  of  the  insulator  by  the  high  field  established  by  the 

0 

positive  charge  deposited  oy  the  impacting  positive  gas  ions. 
Identification  of  the  mechanism  involved  here  may  be  done  by 
use  of  a-c  conductance  measurements  made  as  a  function  of  gate 
voltage  on  the  MOS  structure,^*  and  we  plan  to  pursue  this  line 
of  attach  at  a  later  point  in  our  investigations. 

Breakdown  measurements  made  on  SiO^  films  by  charging 
them  from  a  negative  corona  in  a  nitrogen  atmosphere  indicated 
a  breakdown  field  that  was  approximately  30 %  smaller  than  for 
the  corresponding  charging  in  air,  whereas  for  charging  trom  a 
positive  corona  the  breakdown  fields  were  the  same  for  the  two 
gases.  This  is  a  preliminary  indication  that  for  positive 
surface  potentials  the  current  through  the  silicon  dioxide  is 
controlled  by  the  properties  of  the  silicon/silicon-dioxide 
interface,  presumably  by  the  tunneling  of  electrons  from  the 
silicon  into  the  conduction  band  of  the  silicon  dioxide  where 
they  are  transported  to  the  exposed  surface  without  strong 
trapping.  In  these  experiments,  if  holes  were  injected  into 
the  valence  band  of  the  oxide  (and  this  depends  on  the  impinging 
ion  species,  as  discussed  above)  they  were  rapidly  trapped  at 
or  near  the  exposed  surface.  The  identification  of  the  impinging 
ion  species  in  these  early  experiments  is  not  certain  because 
of  the  important  role  of  impurities  in  the  corona  charging  of 
insulator  surfaces.  Thus  Shahin^*"*  found  that  the  ion  species 
critical  for  surface  charging  by  corona  in  air  are  those 
deriving  from  trace  constituents,  namely  C0^  for  negative 
charging  and  (HgCO^H  positive  charging. 

We  have  concluded  that  the  corona  charging  technique 
has  a  very  high  potential  indeed  for  the  elucidation  of  the 
electronic  behavior  of  the  oxide,  and  other  insulating,  firms 
on  silicon.  However  careful  control  of  the  discharge  conditions 
is  necessary.  Therefore  we  have  had  constructed  an  environmental 
chamber  with  an  air  lock  for  the  introduction  and  removal  of 
samples  and  with  provision  for  controlled  circulation  of  the 
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ambient  gas  through  a  desiccator.  Manipulation  of  tne  sample 
and  of  prcoes  can  be  achieved  from  the  exterior  via  the  use  of 
suitable  gloves  provided  with  the  equipment.  This  chamber 
will  provide  the  required  control  of  the  ambient  gas  and  will 
enable  us  to  mike  a  comprehensive  study  of  eff acts  related 
to  the  propert  es  of  different  species  of  gas  ions.  In 
addition  to  air  and  nitre  gen,  corona  charging  studies  will  be 
made  in  CO^ ,  SFg,  argon  and  neon.  These  studies  will  be  made 
cn  AlgO^  and  Si^  films  on  Si  as  well  as  on  SiOg  films,  and 
r n  composites  of  these  films 

3 •  Locally  Destructive,  Self-Quenching  Breakdown  of  the  Films 

(D.  Y.  Yang  collaborating) 

In  the  self-quenching  technique  of  observing  events 

associated  with  insulator  breakdown,  a  very  thin  metallization 

of  the  top  surface  is  used.  A  local  breakdown  event  in  the 

insulator  causes  the  metal  in  the  immediate  vicinity  of  the 

fault  to  be  removed  by  either  vaporization,  melting  flow  or 

blow-out,  thus  achieving  two  results:  the  damage  to  the 

insulator  is  confined  to  a  very  small  area  (of  micron 

dimension),  and  the  removal  of  the  metallic  electrode  fAom  the 

vicinity  of  the  fault  establishes  an  open  circuit,  which  permits 

one  to  continue  testing  the  remainder  of  the  sampie , ^ *12 

We  have  embarked  on  a  program  of  experiments  to  utilize  this 

technique  systematically  in  the  study  of  breakdown  in  SiO^, 

AI2O3,  and  Si^N^  films  and  their  composites.  We  have  designed 

and  constructed  a  ramp-voltage  generator  which  can  be  cycled 

(i„e.,  produce  breakdowns)  at  rates  in  the  range  of  1  to  100 

per  second.  A  sample  holder  has  been  built  which  provides 

sample  temperatures  in  the  range  of  77  to  500°K 

Preliminary  measurements  have  been  made  at  room 

temperature  on  A1  0  films  pyrolytically  deposited  to  a 
*-  ->  o 

thickness  of  approximately  1000A  on  silicon  substrates  with 
resistivities  of  approximately  1  ohn-cm  and  with  gold  contacts. 
With  the  contact  biased  negatively,  locally  destructive  breakdown 
events  were  observed  at  field  strengths  of  approximately 
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1.1  x  10  V/ cm  for  n-type  substrates  and  8.8  x  10  V/ cm  for 

p-type  substrates.  When  the  gold  contact  was  biased  positively, 

however,  an  unusual  difference  was  observed  between  the  results 

obtained  with  the  two  substrates:  for  an  n-Si  substrate, 

locally  destructive  breakdown  events  were  observed  at  a  field 

strength  of  8.7  x  10^  V/cm,  whereas  for  the  p-Si  substrate 

there  were  no  localized  breakdown  events  at  all,  but,  instead, 

at  field  strengths  approaching  8  x  10  V/cm  the  sample  began  to 

conduct  a  comparatively  large  current  which  thereupon  evaporated 

the  entire  gold  electrode.  Measurement  of  the  MOS  flatband 

voltage  showed,  for  both  n  and  p  substrates,  negative  insulator 

charging  after  breakdown  produced  with  the  geld  contact,  positive, 

and  positive  insulator  charging  after  breakdown  with  the  goid 

contact  negative.  These  results  suggest  that  in  both 

polarities  the  breakdown  is  initiated  by  carrier  injection  from 

the  Si  substrate  into  the  oxide  film. 

The  surface  of  each  sample  was  observed  through  an 
optical  microscope  during  a  sequence  of  breakdowns  Locally 
destructive  breakdown  events  were  observed  to  produce  removal 
of  metallic  electrode  over  areas  with  diameters  of  a  few 
microns.  We  plan  to  examine  samples  which  have  undergone  self¬ 
quenching  breakdown  experiments  in  the  chamber  of  the  Princeton 
high-resolution  scanning  electron  microscope.  This  will  allow 
monitoring  of  the  overall,  as  well  as  the  local,  state  of  the 
sample  as  the  number  of  breakdowns  increase:  the  topography 
of  the  sample,  the  induced  current  image,  the  surface  composition 
(Auger  spectra),  and  the  bulk  composition  (X-ray  fluorescence), 
all  of  which  can  be  recorded  on  film  to  a  high  degree  of 
spatial  resolution. 

The  difficulty  of  achieving  direct  experimental 
observation  of  a  local  breakdown  event  is  immediately  evident 
if  we  note  that:  i)  the  entire  event  is  completed  on  a  time 
scale  of  a  hundred  nanoseconds  or  tens  of  nanoseconds,  or  even 
less,  depending  on  the  particular  structure,  and  ii)  the  lateral 
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spatial  extent  of  the  event  is  on  the  order  of  microns  at  most, 
and  possibly  much  smaller  in  its  earlier  stages,  and  iii)  until 
possibly  well  along  in  its  individual  history,  the  entire 
breakdown  event  is  probably  'buried1  underneath  the  electrode. 

With  regard  to  the  first  point,  sufficiently  fast  electronics 
is  available  to  study  the  time  evolution  of  a  single  breakdown 
event.  Indeed  in  ' constant-voltage '  mods  studies  of  self¬ 
quenching  breakdown  of  SiC>2  films  on  Si,  electrode  potentia1  vs. 

time  (V(t)  vs.  t)  was  followed  during  a  single  breakdown  event 
IS 

by  Klein.  We  believe  that,  by  careful  design  of  the  experiment, 
it  will  be  possible  to  obtain  both  V(t)  vs.  t  and  l(t)  vs.  t 
during  a  single  breakdown  event,  hence  also  I  vs.  V.  Such 
information,  previously  unavailable,  will  be  a  powerful  aid  in 
both  developing  and  checking  theoretical  models  for  the 
breakdown.  Another  line  of  experiments,  also  not  previously 
carried  out,  which  we  expect  to  yield  valuable  information  on 
the  second  point  above  (the  spatial  extent  of  the  local  break¬ 
down  current)  involves  study  of  the  local  breakdowns  in  the 
scanning  electron  microscope  for  a  series  of  almost  identical 
samples,  differing  only  in  the  thickness  of  metallization.  A 
thickness  range  from  approximately  100  to  2000A  will  be  employed. 

A  third  series  of  new  experiments,  also  on  our  agenda  for 
later  undertaking,  involves  the  study  of  laser-induced 
breakdown  in  the  film- 

4.  Charge-Discharge  Studies  of  Trapping  in  the  Films  (N.  Johnson 

collaborating) 

A.  Background  Discussion  and  Experiments 

In  this  portion  of  cur  research  effort  we  are  studying 
the  properties  of  deep  traps  in  the  insulating  films  by  the 
technique  of  first  charging  the  traps  with  carriers  and  then 
selectively  detrapping  the  carriers,  either  by  photon  energy  or 
thermally,  from  traps  that  are  located  progressively  deeper 
in  the  forbidden  gap. 

A  particularly  useful  version  of  the  charge-discharge 
technique,  known  as  the  method  of  Transient  Polarization  Currents, 
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was  originally  used  by  von  Hippel  et  ai  for  the  alkali  halides 
and  later  by  Tabak  and  Warterx^  for  amprphous  selenium.  This 
technique  eliminates  the  need  for  injecting  contacts  and 
requires,  in  fact,  a  blocking  contact,  the  latter  being  far 
easier  to  obtain  on  an  insulator-  The  traps  in  the  insulator 
are  first  charged,  for  example  by  irradiation  of  the  insulator 
with  supra-bandgap  UV  photons  to  create  hole-electron  pairs 
with  subsequent  trapping,  or  alterna^ ively  by  injection  of 
charge  carriers  into  the  insulatoi  by  internal  photoemission 
from  a  contact  or  by  avalanche  injection  from  the  semiconductor 
substrate.  Subsequent  to  the  charging  of  the  traps,  a  step  of 
voltage  is  applied  across  the  insulator.  This  results  in  a 
transient  current  as  the  traps  are  emptied  of  carriers  (it  is 
here  that  the  blocking  contact  is  required).  If  the  voltage 
is  sufficiently  large  so  as  to  sweep  all  free  carriers 
continuously  from  the  sample,  the  resulting  current  for  a 
single  discrete  trapping  level  is  given  by 

qALN 

i(t )  =  -^7 —  e  ,  (1) 

where  q  is  the  magnitude  of  the  electronic  charge,  A  is  the 

area  of  the  blocking  contact ,  L  is  the  thickness  of  the  sample , 

N  is  the  density  of  the  traps,  and  t  is  the  average  time  for 
"t 

emission  of  a  trapped  carrier.  For  purely  thermal  emission 
from  the  traps, 


1 

<ov>  N 


AE/kT 


where  a  is  the  capture  cross  section  of  the  trap,  v  is  the 

carrier  velocity,  <av>  is  the  average  over  the  free-carrier 

distribution  function  of  the  product  av,  N  is  the  effective 

s 

density  of  states  in  the  conduction  band  (assuming  that  the 

emitted  carriers  are  electrons),  AE  is  the  energetic  depth  of 

the  trap  below  the  edge  of  the  conduction  band,  k  is  the 

l6 

Boltzmann  constant,  and  T  is  the  absolute  temperature. 


The  foregoing  expressions  have  been  written  for  a 
single  trapping  level.  In  general.,  there  may  be  a  set  of 
discrete  levels  or  perhaps  a  continuum  of  levels.  To  examine 
these,  the  experiment  is  repeated  at  each  of  a  set  of  successively 
higher  temperatures  to  explore  the  traps  that  lie  progressively 
deeper  in  the  forbidden  gap. 

If  the  sample  is  irradiated  with  light  having  hv  ^  AE 
during  the  progress  of  the  transient,  the  emission  of  carriers 
will  be  enhanced  above  the  thermal  rate,  i.e.,  the  depth  of  the 
trap  is  effectively  reduced  in  value.  This  is  a  particularly 
convenient  technique  for  investigating  traps  which  are 
relatively  deep  in  the  gap,  and  avoids  the  use  of  temperatures 
at  which  ionic  currents  would  become  appreciable.  In  studying 
wide  bandgap  insulators  there  is  no  problem  whatever  in  obtaining 
light  of  sufficient  intensity  to  dominate  the  thermal  emission 
process.  The  minimum  value  of  photon  energy  required  for 
emission  provides  a  further  measure  of  AE.  The  use  of  light, 
instead  of  temperature ,  to  de-trap  the  elect  rons  makes  it 
possible  to  do  the  entire  experiment  at  a  fixed  temperature, 
say  room  temperature  or  even  lower.  The  range  of  trap  energies 
in  the  forbidden  gap  which  can  be  scanned  by  this  optical 
technique  is  far  greater  than  that  which  can  be  scanned  by  the 
conventional  temperature  technique. 

B.  Experimental  Results 

Typical  results  obtained  by  the  use  of  this  technique 
are  shown  in  Figs.  3  and.  4.  Figure  3  shows  the  temperature 
controlled  transient  current  obtained  with  a  4000A  sample  of 
SiOg  on  a  substrate  of  n-Si  with  an  A1  field  plate.  The  traps 
were  first  populated  by  irradiating  the  insulator  with  soft 
(5  kV)  X-rays.  In  this  example,  the  sample  was  held  at  a 
temperature  of  T5°C  and  a  step  of  10  volts  (substrate  negative) 
was  applied.  As  evidenced  by  the  straight-line  segment  on  the 
semi-log  paper,  the  transient  was  exponential  over  most  of  the 
time  range.  (Shallower  traps  can  be  seen  emptying  for  t  <  2  min.) 
Application  of  Eq.  (l)  yields  the  following  results  for  this  se^ 
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of  traps:  Nt  =  1.1  x  10l4  cm-3  and  t  -  2,2  x  10  sec.  If  we 
assume  the  reasonable  values  a  v,  10  cm  ,  v  ^  10  cm/sec, 
and  H  'v  lO1^  cm-3,  the  use  of  Eq.  (2)  yields  the  depth  of  the 
trap:  AE  ^  0.9  eV.  (The  value  of  AE  is  not  particularly 
sensitive  to  the  exact  values  chosen  for  a,  v,  and  Nq.)  The 
value  of  AE  can  be  confirmed  in  two  ways:  either  optically, 
or  from  the  slope  of  log  t  vs.  1/T.  The  confirmed  value  of 
AE  can  be  used  to  check  the  value  used  for  a. 

Another  series  of  experiments  were  performed  on  a 
3500A  sample  of  Si02  thermally  grown  on  an  n-Si  substrate  and 
implanted  with  sodium  ions  to  a  dose  of  10  cm  .  The 
implantation  energy  was  20  kV,  which  produces  a  gaussian 
distribution  of  sodium  atoms  with  the  peak  located  about 
400A  deep  and  with  a  half  width  of  approximately  500A.  In 
this  experiment  the  traps  were  populated  by  the  internal 
photoemission  of  electrons  from  the  substrate  into  the  insulator, 
using  light  of  the  proper  photon  energy  and  a  semitransparent 
field  plate.  The  transient  polarization  current  was  produced 
at  room  temperature,  using  a  bias  of  2k  volts  (substrate 
negative).  The  sample  was  irradiated  with  light  during  the 
progress  of  the  transient,  and  the  intensity  of  the  light 
was  such  that  photon-stimulated  emission  dominated  thermally- 
stimulated  emission.  Determination  of  the  flat-band  voltage 
from  high-frequency  (l  MHz)  C-V  tests  provided  a  measure  of 
the  charge  stored  in  the  insulator.  The  results  of  C-V 
measurements  are  shown  in  Fig.  k.  Curve  1  was  taken  on  the 
sample  as  received.  Curve  2,  which  was  taken  after  charging  the 
traps  with  electrons,  shows  a  positive  flat-band  shift  of 
approximately  6  volts,  indicating  the  deep  trapping  of 
electrons.  In  producing  the  transient  current  using  photon- 
enhanced  emission,  it  was  found  that  a  minimum  photon  energy  of 
approximately  3.5  eV  was  required.  Curve  3  was  obtained  after 
the  transient,  during  which  the  sample  was  illuminated  with 
light  with  the  aforementioned  photon  energy.  The  negative 
shift  of  flat-band  voltage  of  approximately  1.3  volts  indicates 


that  a  fraction  cf  the  traps  released  their  electrons  during 
the  transient  and  are  therefore  located  at  a  depth  of  3.5  eV 
or  slightly  shallower,  but  the  majority  of  the  traps  retained 
their  charge  and  are  therefore  deeper  in  the  gap  than  3-5  eV . 
After  the  current  decayed  to  the  background  level,  shutting 
off  the  ill’imination  and  then  reapplying  it  produced  no 
further  effect;  hence  the  traps  at  the  3.5  eV  depth  must  be 
considered  iv  be  empty.  The  experiment  is  repeatable  in  that 
the  traps  can  be  charged  and  discharged  repeatedly  with 
reproducible  results 

If,  for  simplicity  of  analysis,  we  assume  a 

rectangular  distribution  of  traps  that  extend  to  a  depth  of 
o 

600A  into  the  oxide  (corresponding  roughly  to  the  distribution 

of  sodium  atoms),  computation  shows  that  the  flat-band  shift 

from  Curve  2  to  Curve  3  of  Fig-  4  indicate^  a  density  of  traps 

of  approximately  10  cm  at  a  depth  of  _j.5  eV  or  slightly 

shallower.  This  corresponds  to  a  number  of  traps  per  unit  area 

12  -2 

that  is  slightly  less  than  10  cm  .  This  is  considerably 

15  -2 

smaller  thai.  the  dosage  of  sodium  ions,  which  was  10  cm 

At  this  time  it  is  a  matter  of  conjecture  whether  the  electron 

traps  observed  here  are  due  to  the  displacement  damage  caused 

by  implantation  or  are  caused  by  the  sodium  atoms  directly.  If 

it  is  the  latter,  only  a  small  fraction  of  the  sodium  atoms 

are  active  in  the  trapping  process  Note  that  with  an 

15  2 

implantation  dose  of  10  '  Na  lons/cm  a  maximum  of  approximately 

TO  p 

1% ,  namely  10  J/cm  ,  can  remain  unneutralized,  even  during  the 

implantation.  This  is  because  the  oxide  breakdown  field 

(5  x  106  to  10?  V/cm)  is  reached  at  the  latter  concentration. 

Thus,  field -induced  electron  injection  into  the  oxide  must  take 

place  throughout  virtually  the  entire  period  of  the  implantation 

This  is  a  complication  requiring  further  investigation. 

Examples  of  results  obtained  using  another  version 

of  the  charge-discharge  technique  are  shown  in  Figs.  5  through  7 

The  sample  here  consisted  of  a  film  of  SiO  thermally  grown  to 

o  e- 

a  thickness  of  approximately  1000A  on  an  n-Si  substrate  and 
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subsequently  implanted  with  aluminum  ions  to  a  dose  of  10li;  cm-2. 

The  implantation  energy  was  20  keV,  which  produces  a  gaussian 

distribution  cf  aluminum  atoms  with  the  peak  located  approximately 

at  the  middle  of  the  insulating  film  and  with  a  half  width  of 
o 

about  500A.  In  Fig  5  is  shown  a  series  of  high-frequency 

(l  MHz)  C-V  curves.  Curve  1  is  for  the  sample  as  received.  The 

flat -band  voltage  of  -15  volts  indicates  an  initial  storage 

of  positive  charge  in  the  insulator  equivalent  to  3  x  1012 
2 

cnarges/ cm  if  the  charge  is  located  at  the  silicon-insulator 
interface,  and  double  this  amount  if  the  center  of  gravity  of 
the  charge  is  located  at  the  middle  of  the  insulator.  The 
presence  of  deep  electron  traps  in  the  insulator  is  shown  by 
Curves  2  through  7  These  curves  were  produced  in  succession. 

For  each,  the  gate  voltage  was  held  at  a  constant  positive 
value  and  electrons  were  photoin j ected  from  the  silicon  into 
the  insulator  where  they  were  deeply  trapped.  It  will  be  seen 
that  the  flat-band  voltage  of  the  MOS  structure  is  very  nearly 
equal  to  the  voltage  at  which  the  gate  was  held  during  the 
photoinjection  of  electrons.  Presumably  the  trapping  of 
electrons  proceeds  until  the  negative  charge  in  the  insulator 
is  sufficient  to  reduce  the  electric  field  at  the  silicon- 
insulator  interface  to  zero,  thus  bringing  to  a  halt  the 
injection  of  further  electrons  Furthermore,  the  electron 
traps  do  not  saturate  with  gate  voltages  of  up  to  35  volts. 

This  puts  a  lower  limit  on  the  number  of  traps:  1.5  x  1013  cm-3 
if  the  center  of  gravity  of  the  charge  is  approximately  at  the 
middle  cf  the  insulator.  Figure  6  shows  the  transient  current 
that  flowed  in  the  trap-charging  process  that  took  the  C-V 
characteristics  of  the  structure  from  Curve  1  to  Curve  2  in 
Fig.  5.  Figure  7  shows  the  recovery  l approximately  exponential 
in  time)  of  the  flat-band  voltage  of  the  sample  when  it  is 
irradiated  with  light  from  a  xenon  source  with  the  field 
plate  of  the  sample  biased  negatively-  Further  experiments  are 
underway  to  determine  more  accurately  the  density  of  the  traps 


and  to  measure  their  capture  cross  sections  and  depth  in  the 
forbidden  gap.  The  results  obtained  to-date  show  the  effective¬ 
ness  of  the  charge-discharge  techniques  for  obtaining  valuable 
information  on  the  electronic  behavior  of  thin  oxide  films. 

A  brief  summary  of  these  results  is  being  prepared  for  publi¬ 
cation. 

5-  Electron-Beam  Probe  Studies  of  the  Films  (C,  T.  Shih 
collaborating) 

In  this  set  of  experiments  ,  the  problem  of  creating 
a  suitable,  equivalently  ohmic  contact  to  the  insulator  is 
dealt  with  by  using  a  beam  of  electrons  to  create  a  reservoir 
of  electrons  just  under  the  surface  of  the  insulator.  A 
front  electrode,  of  very  thin  metal  and  comparatively  transparent 
to  the  electron  beam,  is  used  to  impress  voltage  across  the 
sample.  The  energy  of  the  electron  beam  can  be  adjusted  in  the 
range  of  1-5  kV. 

The  operation  of  the  experiment  can  be  visualized  with 
the  aid  of  the  upper  band  diagram  shown  in  Fig,  8.  The  incident 
high-energy  electrons,  after  penetration  of  the  thin  field  plate, 
generate  hole-electron  pairs  in  the  insulator  up  to  their  depth 
of  penetration  R.  (Also  produced  will  be  X-rays  which  can 
penetrate  the  entire  sample.)  In  SiO^  the  holes  have  a  very 
small  mobility  and  are  trapped  and  ultimately  recombined  with 
electrons  near  their  point  of  formation.  In  the  vicinity  of 
the  cathode,  the  density  of  electrons  exceeds  the  density  of 
holes  because  of  the  electrons  injected  into  the  insulator 
by  the  primary  (incident)  beam  directly  as  well  as  by  a 
secondary  current  generated  by  the  primary  beam  in  the  gate  by 
collisional  excitation.  The  resulting  negative  space  charge 
creates  a  maximum  in  electron  potential  energy  as  shown  in  the 
diagram.  (This  causes  the  effective  sample  voltage  to  be 
larger  than  the  applied  voltage,  V  ,  and  also  causes  the  effective 
sample  length  to  be  smaller  than  the  actual  length,  as  is 
indicated  in  the  diagram.  This  effect  must  be  taken  into  account 


Fig.gi,  High-field  Separation  of  £Uclron'ho/e  pairs 
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in  the  analysis  of  the  data).  In  Si02  electrons  are  generally 
not  strongly  trapped  and  therefore  application  of  the  voltage 
across  the  sample  results  in  their  transport  through  the 
insulator  to  the  substrate,  which  serves  as  the  anode.  Figure  9 
shows  the  sample  current  as  a  function  of  sample  voltage  for 

O 

Sample  HH71*,  which  consisted  of  U000A  of  Si02  thermally  grown 
on  a  substrate  of  n-Si .  Here  the  beam  voltage  was  1.6  kV  and 

7 

the  beam  current  was  1.0  x  10  amp  It  can  be  seen  that  for 

sample  currents  considerably  smaller  than  the  beam  current, 

the  sample  current  follows  approximately  a  sqaure  law.  This 

is  followed  by  a  steeply  rising  section  beyond  which  the 

current  gain  exceeds  unity ,  that  is ,  the  sample  current 

exceeds  the  beam  current.  A  possible  explanation  for  this 

effect  is  indicated  by  the  lower  band  diagram  of  Fig.  8.  Here 

the  voltage  applied  to  the  sample  has  overcome  the  effect  of 

the  negative  space  charge  in  the  region  of  the  cathode,  and  the 

holes  that  have  oeen  created  less  than  one  Schubweg  from  the 

contact  (0  <  x  <  L  in  Fig.  8)  are  swept  out  and  augment  the 
s 

sample  current.  Inasmuch  as  one  high-energy  electron  from  the 
beam  will  create  many  hole-electron  pairs  in  the  insulator, 
the  sample  current  can  exceed  the  beam  current  if  holes  are  able 
to  leave  via  the  contact  instead  of  recombining  within  the 
insulator.  If  this  explanation  for  the  current  gains  exceeding 
unity  were  the  only  possible  one  it  would  be  of  considerable 
importance,  since  it  would  directly  establish  hole  motion  in 
the  oxide  and  provide  detailed  quantitative  information  about 
it.  Actually,  it  is  highly  likely  that  some  of  the  current 
gain  has  its  origin  in  the  metal  gate  as  already  mentioned 
above.  The  thickness  of  metallization  employed  to-date  is  such 
that  a  substantial  fraction,  up  to  one-half,  of  the  incident 
beam  energy  is  lost  in  the  metal  gate  producing  a  large  number 
of  hot  electrons  in  the  gate.  Some  of  the  latter  will  escape 
into  the  oxide  providing  a  secondary  contribution  to  the  total 
electron  current  delivered  to  the  oxide.  In  order  to  derive 
any  conclusions  about  hole  motion  in  the  oxide  it  is  necessary 
to  eliminate  or  separate  out  the  secondary-current  contribution 
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of  the  gate .  Towards  this  end  we  plan  to  prepare  a  series  of 

o 

samples  with  ultra-thin  metallizations,  down  to  100A  thickness. 

Further  important  results  are  shown  in  Figs.  10  and 

11.  Figure  10  shows  the  dependence  of  the  sample  current  on  the 

beam  cuirent  in  the  square-l£.w  regime  of  Is  vs.  Vs  ,  and  indicates 

that  the  sample  current  is  approximately  linearly  dependent  on 

the  magnitude  of  the  beam  current  (the  lack  of  parallelism 

among  the  curves  is  attributable  to  the  difference  between  the 

applied  and  effective  sample  voltages;  see  Fig.  8).  The  results 

shown  in  Fig.  11  indicate  that  the  sample  current  is  only 

weakly  dependent  on  the  temperature.  Finally,  analysis  of  the 

data  in  the  square-law  regime  shows  that  the  sample  current  is 

considerably  smaller  than  the  magnitude  predicted  by  the  Mott- 

Gurney  relationship  for  single-carrier  space-charge-limited 

2  3 

currents  in  a  trap-free  insulator:  J  =  (9/8)eyVs/L  ,  where  e 

is  the  absolute  dielectric  constant  of  the  insulator  and  y  is 

the  mobility  of  the  charge  carrier.  Taken  together,  these  data 

lend  themselves  to  the  following  interpretation:  The  current 

is  predominantly  single-carrier  and  space-charge-limited,  but 

a  large  fraction  of  the  electrons  are  trapped  and  the  appropriate 

IT 

expression  for  the  current  is,  consequently, 

V2 

J  =  |  0ey  ,  (3) 

8  I,3 

where  0  is  the  ratio  of  free  to  trapped  electrons,  and  0  <<  1. 

The  comparative  independence  of  current  on  temperature  and  the 
nearly  linear  dependence  on  magnitude  of  beam  current  indicates 
that  the  emission  of  charge  carriers  from  traps  is  not  thermally 
produced  but  is  instead  the  result  of  a  photon  flux  which  is 
generated  over  the  range  of  the  incident  primary  beam,  that 
is,  in  the  gate  and  in  the  first  few  hundred  Angstroms  of  the 
oxide  near  the  gate. 

In  order  to  investigate  the  foregoing  effects 
further,  we  are  currently  arranging  the  equipment  for  illumination 
of  the  sample  with  light  from  a  high-intensity  monochromator. 
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Thi.  will  permit  us  tc  study  the  effect  on  6  of  the  intensity 
and  photon  energy  of  the  light  and  will  provide  a  method  of 

investigating  the  density  of  traps  and  their  distribution  in 
energy. 

6'  ^g£onflary-Electron  Spectrometry  (D.  Guterman  collaborating) 
A-  Background  Discussion 

ie  pas  few  years  secondary  electron  spectroscopy 
has  become  a  powerful  tool  for  the  analysis  of  solid  surfaces.  ' 
TT'WC  ~eParatg  analytical  techniques  have  been  developed  for  this 
purpose  using  a  spectroscopy  of  secondary  electrons.  The  first 
is  Auger  electron  spectroscopy  (AES)  designed  to  determine  the 
elemental  composition  of  a  surface.  The  second  is  electron 
spectroscopy  for  chemical  analysis  (ESCA)  which  is  designed  to 
determine  the  chemical  environment  in  which  the  atomic  con- 
-Utuents  of  a  surface  reside  The  construction  of  a  secondary 
electron  spectrometer  with  both  high  resolution  and  high 
ensitivity  will  allow  an  extension  of  AES  to  perform  analysis 
of  the  chemical  environment  of  the  atoms  on  a  solid  surface 
and  to  probe  the  band  structure  near  the  surface,  much  as  is 
done  in  external  photoemission  studies.  In  AES  the  sample  is 
bombarded  with  an  electron  beam  of  a  few  thousand  electron 
voxts  energy  and  the  distribution  of  secondary  electrons  is 
measured.  The  energies  at  which  peaks  occur  in  the  secondary 
emission  spectra  correspond  to  energy  level  differences  in  the 
emitting  atom.  Therefore  the  elemental  species  may  be  identified 
through  a  measurement  of  the  peak  energies.  The  amplitude  of 
these  peaks  above  the  background  level  is  a  measure  of  the 
relative  abundance  of  the  various  atomic  species  found  on  the 
surface-  This  technique  has  a  sensitivity  in  the  range  of  10"2 
tC  bet+er  than  10~3  monolayers  depending  on  the  shot  noise  in 
the  particular  detector  system,  the  particular  elements  being 
measured, and  the  elemental  composition  of  the  background  matrix. 

In  addition  to  this  powerful  capability  of  identifying  atomic 
species  on  a  surface  there  is  a  great  deal  of  information 
stored  In  the  fine  structure  of  the  Auger  peaks  which  has  yet  to 
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be  fully  explored  and  utilized.  The  energy  of  an  Auger  pea* 
is  determined  by  the  energy  level  differences  within  the  atom 
and,  particularly  where  valence  orbitals  are  involved,  these 
differer-  is  are  altered  by  the  chemical  bonding  configuration 
of  an  atom.  Therefore  small  shifts  in  energy  of  a  particular 
peak  can  be  related  to  the  chemical  bonding  conditions  of  the 
atom  being  probed.  The  fine  structure  on  the  low  energy  side 
of  the  Auger  peaks  and  the  primary  electron  beam  can  be 
associated  with  plasmon  and  phonon  excitations.  A  spectrometer 
with  sufficient  energy  resolution  would  therefore  also  have  the 
capability  of  measuring  the  plasmon  and  phonon  spectra  near  the 

surface. 

The  shape  of  the  Auger  peaks  is  determined  in  large  part 
hy  the  density  of  states  in  the  valence  hand  near  the  surface. 

This  is  due  to  the  fact  that  the  probability  of  an  Auger  electron 
of  energy  E,  to  he  emitted  in  a  WXY  transition  (W  denotes  the 
hand  where  the  initial  vacancy  is  formed  hy  the  incident  electron 
beam,  X  the  hand  of  origin  for  the  electron  decaying  to  fill  the 
W  hand  vacancy,  and  Y  the  hand  of  origin  for  the  electron 
ejected  to  relax  the  energy  difference  between  hands  W  and  X) 
is  proportional  to  the  density  of  filled  states  in  the  X 
and  Y  hands.  Experimentally  this  means  that,  knowing  the 
distribution  of  electrons  in  an  Auger  peak,  one  can  determine 
the  density  of  filled  states  in  the  hands  participating  in  the 
emission  process.  For  transitions  of  the  type  KW  or  LW 
where  the  valence  bund  is  involved  we  may  determine  the  electron 
density  of  states  in  the  valence  hand.  The  analysis  is  simplified 
in  these  two  cases  hy  the  fact  that  the  K  and  L  hands  are 
usually  very  narrow  and  therefore  may  he  considered  mono-energetic. 

The  benefits  to  he  derived  from  an  appropriately  designed 
secondary  electron  spectrometer  and  a  highly  focused  electron 
beam  for  exciting  the  secondaries  may  he  summarized  as  follows: 

1.  Elemental  analysis  of  a  surface  may  he  carried  out  with 
a  sensitivity  of  better  than  iCf3  monolayers  under  optimum  conditions. 
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2-  The  chemical  bonding  state  of  the  elements  on  the  surface 
of  a  solid  may  be  determined  with  essentially  the  same  sensitivity 
as  the  elemental  analysis. 

3  A  measurement  of  the  density  of  occupied  states  near 
the  surface,  for  bands  participating  in  the  Auger  emission,  may 
be  carried  out. 

h.  High  spatial  resolution  for  the  analytical  capabilities 
in  1-3  above  may  be  achieved.  The  resolution  is  limited  only 
by  the  size  of  the  incident  electron  beam. 

B.  Secondary-Electron  Spectrometer 

In  order  to  obtain  all  of  these  benefits  in  an  electron 
spectrometer  we  have  designed  a  system  which  has  the  maximum 
possible  collection  efficiency  and  throughput  for  the  secondary 
electrons  and  the  lowest  possible  shot  noise.  In  addition,  the 
spectrometer  was  designed  with  a  high  energy  resolution  so  that 
the  measurement  of  fine  structure  in  the  Auger  peaks  as  well 
as  an  accurate  measurement  of  the  peak  positions  would  be 
possible  The  system  resulting  from  these  design  considerations 
is  shown  in  Fig.  12.  High  collection  efficiency  is  accomplished 
since  the  spectrometer  entrance  subtends  an  angle  of  96  degrees 
with  respect  to  the  sample  The  throughput  is  limited  only  by 
the  transparency  of  the  first  two  grids  and  is  maintained  at  a 
high  level  by  using  a  wire  mesh  for  the  grids  with  an  optical 
transparency  of  90  per  cent.  The  functions  of  individual 
components  of  the  spectrometer  are  as  follows: 

1.  Grid  one  is  grounded  during  normal  operation  to  provide 
a  field  free  region  around  the  sample  and  thereby  prevent  errors 
in  the  Auger  energies  caused  by  fields  between  the  sample  urd 
the  detector  entrance.  To  increase  the  collection  of  secondary 
electrons,  when  high  resolution  is  not  necessary,  grid  one  may 
be  biased  positive  with  respect  to  the  sample. 

2.  Grid  two  is  biased  negatively  with  respect  to  the  sample 
by  a  ramp  power  supply.  This  provides  a  low  energy  cut-off 
filtering  of  the  secondary  electrons.  The  wire  mesh  used  to 
fabricate  this  grid  is  very  fine,  resulting  in  a  theoretical 


Fig.  12  Diagram  of  the  Secondary  Electron  Spectrometer 
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resolution  of  0,01  percent  of  the  energy  of  the  secondaries. 

3.  Grid  three  is  held  at  a  potential  slightly  more  negative 
than  that  of  grid  two.  This  has  two  beneficial  effects;  the 
first  is  the  focusing  of  low  energy  electrons  passing  through 
grid  two  toward  the  aperture  in  grid  four.  The  second  is  to 
prevent  secondaries,  generated  on  the  walls  of  the  detector  body 
by  high  energy  electrons  reflected  from  the  sample,  from 
entering  the  spectrometer,  thereby  keeping  the  shot  noise  to  a 
minimum. 

h ,  Grid  four  is  held  exactly  at  the  potential  of  grid  two. 
This  provides  a  low  field  region  where  electrons  passing  through 
the  second  grid  are  focused,  by  the  collector  potential  and 
the  potential  on  grid  three,  toward  the  collector.  Electrons 
which  pass  through  the  second  grid  with  high  energy  are  not 
deflected  to  the  collector  by  this  low-strength  focusing  field 
and  therefore  the  necessary  high  cut-off  filter  is  established. 
This  minimizes  the  current  through  the  detector  and  thereby  the 
shot  noise  of  the  spectrometer. 

5.  The  shield  between  the  sample  and  grid  one  prevents 
high-energy,  reflected  electrons  from  reaching  the  collector  and 
thereby  contributing  to  the  shot  noise  of  the  system. 

6.  The  collector  is  a  high  gain,  low  noise,  continuous, 
dyncde  electron  multiplier.  The  potential  at  the  front  of  this 
detector  is  positive  with  respect  to  grid  two  to  draw  low 
energy  electrons  into  the  multiplier  and  it  may  be  varied  to 
provide  a  real-time  control  of  the  resolution-sensitivity  trade¬ 
off.  The  potential  drop  across  the  multip  Ler  is  controlled 
separately  to  ensure  the  highest  possible  gain  which  does  not 
saturate  the  multiplier  under  various  signal  strength  conditions. 

7.  The  outer  body  of  the  detector  provides  structural 
strength  to  the  grid  structure,  prevents  electrons  scattered  in 
the  vacuum  chamber  from  entering  the  spectrometer,  and  provides 
the  magnetic  shielding  essential  to  the  efficient  operation  of 
the  spectrometer. 

The  response  of  this  system  has  been  modeled  on  the 
IBM  360  model  91  computer  and  the  design  has  been  optimized  for 
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resolution  and  collection  efficiency  The  results  of  this 
modeling  indicate  that  the  resolution  under  optimum  conditions 
may  be  as  hign  as  0.025  electron  volts.  In  practice  the 
resolution  is  determined  by  the  second  grid  and  the  shot  noise 
which  is  limited  by  the  high  energy  cut-off.  Figure  13  shows  a 
plot  of  the  trajectory  of  the  electrons  passing  through  the 
spectrometer  and  it  illustrates  the  high  energy  cut-off  properties 
of  the  spectrometer.  Figure  lb  shows  the  calculated  efficiency 
of  the  spectrometer  as  functions  of  the  potential  of  the  third 
grid  and  the  front  of  the  electron  multiplier  with  respect  to 
the  grid-two  potential.  It  illustrates  another  valuable  feature 
of  this  electron  spectrometer,  namely  that  the  resolution- 
sensitivity  trade-off  may  be  controlled  by  altering  the  potentials 
in  the  spectrometer  without  any  requirement  for  alteration  of 
the  physical  arrangement  as  is  necessary  with  the  conventional 

cylindrical  mirror  system. 

C.  Preliminary  Results 

The  spectrometer  and  an  appropriate  power  supply  have 
been  completed.  Preliminary  spectra  have  been  taken  with  this 
spectrometer  and  as  a  standard  they  have  been  compared  with  the 
spectra  from  a  commercial  (Varian  Associates)  LEED-Auger  device. 

At  the  present  time  the  new  system  is  exhibiting  a  signal-to- 
noise  ratio  more  than  fifty  times  better  than  the  commercial 
device  under  identical  conditions  of  excitation.  Figures  15  and 
16  show  typical  spectra  taken  from  the  new  spectrometer  and  the 
LEED-Auger  device  using  the  same  electron  gun  in  the  same 
vacuum  chamber.  Figure  15  shows  a  spectrum  taken  with  the 
LEED-Auger  device  with  250  microamps  incident  beam  current 
and  a  modulation  (necessary  for  differentiation  of  the  signal 
with  a  lock-in  amplifier)  of  10  volts  peak-to-peak.  These 
parameters  represented  the  minimum  for  obtaining  a  usable 
spectrum.  Figure  l6  shows  a  spectrum  taken  with  the  new  detector 
with  30  microamps  incident  beam  current  and  a  modulation  of 
5  volts  peak-to-peak.  These  changes  in  the  parameters  reduced 
the  signal  by  more  than  a  factor  of  six  and  the  signal -to- 
noise  ratio  has  improved  approximately  by  a  factor  of  eight. 
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Auger  spectrum  of  a  Si  surface  taken  with  a  commercial  LEED-Auger 
device.  Incident  current  250  microamps,  10  volts  peak  to  peak 
modulation. 


Figure  16 


Auger  spectrum  of  a  Si  surface  taken  with  the  new  spectrometer 
Incident  current  80  microampsl  5  volts  peak  to  peak  modulation 
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This  system  will  be  applied  to  the  study  of  the  elemental 
distribution,  chemical  composition,  and  surface  band  structure 
of  thin  insulating  films. 

The  spatial  resolution  of  the  electron  spectrometer 
described  above  is  limited  by  the  size  of  the  exciting  electron 
beam.  Therefore,  in  order  to  attain  the  best  possible  spatial 
resolution,  this  system  is  to  be  mounted  on  the  specimen  chamber 
of  a  scanning  electron  microscope  (SEM).  Since  this  technique 
is  essentially  a  surface  analysis,  due  to  the  short  elastic 
range  of  the  secondary  electrons,  it  is  necessary  that  the 
experimental  chamber  provide  an  environment  in  which  the  nature 
of  the  specimen  surfaces  can  be  maintained  constant  for  the 
duration  of  a  measurement.  Commercially  available  SEM  equip¬ 
ment  does  not  have  the  required  vacuum  capability  and  therefore 

=uch  a  system  has  been  built  at  Princeton  University.  At  present, 
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a  system  capable  of  10  7  tcrr  in  the  specimen  chamber  has  been 
completed  and  a  modification  is  presently  underway  to  provide 
this  system  with  an  air  lock  so  that  the  SEM  can  be  maintained 
at  constant  high  vacuum,  except  for  infrequent  maintenance 
requirements.  The  specimen  chamber  is  equipped  with  X-ray 
fluorescence  spectroscopy,  induced  current  probes,  transmitted 
electron  probes,  voltage  contrast  capability,  as  well  as  the 
traditional  detector  for  topographic  imaging.  Measurements 
on  silicon  dioxide  have  been  performed  in  the  SEM  using  the 
Auger  spectrometer  but,  as  yet,  the  system  has  not  been  perfected 
to  allow  the  observation  of  band  structure  in  the  specimen.  An 
example  of  these  preliminary  spectra.  Fig.  17,  shows  that  the 
surface  is  composed  of  silicon,  oxygen  and  carbon.  The  silicon-to- 
oxygen  ratio  was  found  to  differ  slightly  from  one  part  of  the 
sample  to  the  next  although  at  the  present  time  we  have  no 
quantitative  measure  of  this  observation.  The  quality  of  the 
analytical  information  available  should  be  greatly  improved 
by  the  completion  of  the  system  within  the  next  few  months. 

This  SEM  system  has  been  used  to  study  pyrolytic 

o 

A1203  thin  films  v900  to  lUOOA  thickness)  grown  on  p-type  (100) 


silicon  substrates.  Utilizing  the  standard  secondary  electron 
topographic  image,  these  samples  were  found  to  be  smooth  and 
flat  out  to  the  limits  of  detectability  over  most  of  their 

surface.  This  places  an  upper  bound  on  the  topographic  varia- 

o 

ticns  on  the  surface  of  the  order  of  10CA.  There  were  several 
xarge  defects  found  distributed  ever  the  surface  of  each  of 
the  pyrolytic  A^O.^  samples  which  were  of  two  distinct  types. 

The  first  was  a  round  hole  with  pointed  crystallites  rising 
from  the  center  and  the  second  was  a  large  l, up  to  one  micron 
in  diameter),  irregular  shaded  piece  of  crystalline  material. 
There  was  no  spatial  correlation  between  these  defects  or  any 
other  measureable  parameter  on  the  samples.  The  identification 
of  the  material  centered  in  the  round  holes  shown  in  Fig.  l8 
was  accomplished  by  X-ray  fluorescence.  The  material  is 
predominantly  silicon  and  it  is  likely  some  artifact  of  the 
film  fabrication  process  Similarly,  the  irregular  crystalline 
material  found  on  the  same  samples  was  determined  to  be  aluminum, 
also  interpreted  as  an  artifact  of  the  growth  process. 

Additional  Al^O^  samples  are  being  prepared  under  different 
conditions.  Study  of  these  new  samples  should  allow  us  to 
determine  the  origin  of  the  defects. 

D.  Optical  Instrumentation 

In  conjunction  with  the  electron  spectrometry  studies 
of  thin  films,  optical  studies  of  both  internal  and  external 
photoemission  will  be  carried  out  A  high  energy  light  source 
has  been  set  up  and  is  presently  operational.  It  is  composed 
of  a  high  intensity  vacuum  ultra-violet  discharge  lamp  of  the 
Hmteregger  type, a  Specs  1500  vacuum  ultraviolet  high  intensity 
monochromator  (6  8  op+ics),  and  a  LiF  window  and  adapter  to 
various  sample  vacuum  systems  with  a  half  silvered  mirror  in  the 
light  path  for  calibration.  This  light  source  and  monochromator 
combination  provides  m  excess  of  1000  times  the  intensity 
previously  available  at  Princeton  for  the  wavelengths  useful 
in  this  study.  The  light  source  monochromator  system  have  been 
calibrated  in  wavelength,  intensity  and  long  term  stability  and 


V 
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Figure  18 

SEM  micrograph  of  an  Al^O^  film  on  a  silicon  substrate  showing 
holes  in  the  layer  in  which  silicon  crystals  are  growing. 


*5- 


the  output  curve  attached  in  Fig.  18b  shows  high  intensity  out 
to  1^00  angstrom  units.  The  capability  of  this  system  will  be 
applied  to  the  characterization  of  thin  films  of  AlgO^  ,  SiOg  , 
and  Si^N^  in  the  coming  months. 

7 •  Monte  Carlo  Studies  of  Hot-Electron  Behavior  in  Insulating 
Films  at  High  Fields  (S.  Baidyaroy  and  B.  Zee  collaborating) 

A.  Background  Discussion 

In  order  tc  obtain  an  understanding  of  the  physical 
mechanisms  underlying  electrical  breakdown  of  thin  insulating 
films  at  high  applied  fields  it  is  necessary  to  obtain  a  reasonable 
description  of  hct-electrcn  behavior  at  the  high  fields.  This 
is  an  extremely  difficult  theoretical  problem  which,  in  the 
context  cf  the  thin-insulating-film  problem,  is  completely 
beyord  the  reach  of  analysis.  We  have  therefore  sought 
efficient  means  of  utilizing  the  high-speed  digital  computer 
to  study  the  problem.  The  costs  of  computation  have  provided 
a  strong  stimulus  to  find  truly  simple,  namely  one-dimensional, 
models  to  simulate  the  main  features  of  the  full  three- 
dimensional  problem. 

The  problem  we  have  addressed  ourselves  to  is  as  follows: 

An  electron  with  a  specified  initial  energy  is  injected  into  the 

film  in  a  randomly  chosen  direction  at  one  end  of  the  film 

taken  at  z  =  0.  The  electron  is  driven  by  the  field  toward  the 

other  end  of  the  film,  at  z  =  t,  where  t  is  the  film  thickness. 

In  moving  through  the  film  the  electron  picks  up  energy  from 

the  applied  field,  but  it  is  also  randomly  scattered  by 

collisions  with  phonons  (quantized  'lattice'  vibrations  in  the 

film).  It  is  assumed  that,  at  the  high  fields  of  interest  in 

6  7 

the  breakdown  problem,  namely  in  the  range  10  to  10  V/cm, 
the  only  significant  scattering  events  are  optical -phonon 


emissions,  the  fixed  energy 
(An  electron  with  energy  E 


e  being  lost  in  each  such  collision, 
ph 

<  e ^  is  not  permitted  to  make  an 


energy-loss  collision).  The  field-driven  electron  is  followed, 


as  a  classical  particle,  in  its  randomly-colliding  course  through 

18 

the  film,  using  the  Monte  Carlo  (random  number)  method  of 


Fig.  181? 

Output  spectra  of  the  specs  1500  vacuum  UV  monochromator 
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dealing  with  the  probabilistic  features  of  its  trajectory- 
history.  At  the  exiting  end  of  the  film,  at  z  =  t,  the  energy 
and  angle  of  emergence  are  noted  (as  well  as  the  lateral  distance 
traversed  by  the  electron).  A  second  electron  is  then  injected 
and  followed  until  emergence  at  the  plane  z  =  t,  then  a  third 
electron  and  so  on,  up  to  200  or  more  electrons.  In  this 
way  a  histogram  of  emerging  electron  energies  is  complied 
(likewise,  one  for  emerging  angles).  These  particular  histograms 
are  of  special  interest  because  experimental  determination  of 
emergent  electron  energies  is  feasible  -  as  compared  to  the 
extraordinary  difficulties  of  experimental  determination  of 
hot -electron  distributions  inside  a  medium.  A  considerable 
number  of  situations  have  been  studied  on  the  computer 
in  this  way,  mostly  const ant-mean-free-path  (independent  of 
energy)  problems  with  different  scattering  laws,  corresponding 
to  isotropic  scattering,  dominant  (polar  optical  phonon) 
forward  scattering.  In  order  to  simplify  the  computer  calculations, 
in  particular  to  avoid  numerically  solving  an  integral  equation 
over  and  over  again,  we  have  utilized  the  ingenious  self- 
scattering  method  of  Rees,  ^  explained  below.  Further,  we  have 
systematic ally  explored  one— dimensional  simulations  of  the  full 
three-dimensional  problems,  and  have  shown  that  these 
simulations  are  extremely  useful  in  that  they  reproduce  the 
salient  features  of  the  latter.  We  have  further  shown  that 
strong  forward  scattering,  by  itself,  is  insufficient  to  lead 
to  instability  of  the  hot-electron  distribution,  that  is,  to 
'runaway'  of  the  distribution.  Finally,  we  have  initiated  the 
study  of  the  practically  important  class  of  problems  in  which 
the  mean-free -path  is  energy  dependent.  If  the  mean-free-path 
increases  fast  enough  with  energy  we  would  expect  runaway  of 
the  distribution,  and  we  have  shown  that  this  is  indeed  so. 

Further,  we  have  shown  thet  the  problem  of  distribution 
stability  can  be  studied  graphically  using  the  already-obtained 
results  of  the  far  simpler  constant-mean-free-path  computations. 
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B.  Mathematical  Formulation 

The  mean  free  path  (mfp)  X  is  defined  such  that  the 
probability  of  scattering  in  the  increment  of  path  length  ds 
is  ds/X.  For  constant  mfp  (independent  of  electron  energy) 
this  leads  to  the  decay  lav:  n  -  nQ  exp  ( — s/ X )  where  n  is 
the  electron  concentration  (or  number)  after  the  distance 
s  along  the  trajectory  has  been  traversed  and  n^  the  corre¬ 
spond?-  concentration  (or  number)  at  s  =  0,  Defining  P(s) 
as  the  probability  that  the  electron  has  gone  the  distance  s 
without  scatte  '.ng,  we  readily  obtain  the  equation 


dP 

P 


(4) 


which  has  the  solution 


P(s) 


(5a) 


which,  for  constant  mfp  becomes 

P(s)  »  exp(- 


(5b) 


Defining  P  (s)  as  the  probability  density  that  the 
sc 

electron  scatter  in  incremental  path  length  ds  after  traversing 
path  length  s  without  scattering,  we  have 


Pjs)  = 

SC 


^  =  *  exp(:  I.'  r, 


which  for  constant  mfp  reduces  to 


P  (s)  ■  1 
sc  —  exp 


K) 


(6a) 


(6b) 


Note  that  P^s)  is  properly  normalized,  that  is,  that 

1  ds  P  (s)  =  1.  Also  note  that,  for  constant  mfp,  s  P  (s)  ds 

i  sc  sc 

'  c  *  0 

The  correspondence  with  the  random  number  variable  r' 

lying  between  0  and  1  is 


dr'=P(s)ds  ,  (7) 

sc 

which  integrates  to 

rr '  rs 

dr  =  r'  =  P  (s' )  ds'  ,  (8) 

i  I  sc 

"  o  *  o 

Using  a  suitable  computer  routine  to  generate  the  random 
number  r'  ,  the  corresponding  trajectory  length  s  between  successive 
collisions  is  determined  by  solution  of  the  integral  equation  (8). 
In  the  general  case  this  integral  equation  will  be  analytically 
intractable  and  will  have  to  be  solved  numerically  (over  and  over 
again).  A  very  important  exception  is  the  constant  mfp  problem, 
for  which  (8),  using  (6b),  has  the  simple  analytic  solution 

r'  *  1  *  ex5(-  1)  .  (9) 


or,  letting  1-r*  =  r, 

(constant  mfp)  s  =  -X  in  r  ,  (10) 

where  r  is  also  a  random  variable  in  the  interval  0  to  1.  Thus, 
for  the  constant  mfp  problem  s  is  given  simply  by  (10),  once  the 
computer  has  generated  r. 

Without  further  simplification  a  formidable  computer 
calculation  still  remains.  Consider  the  relatively  simple  hot- 
electron  problem  of  constant  mfp  and  isotropic  scattering.  It  is 
easy  to  see  that,  because  of  the  assumed  scattering  isotropy,  this 
problem  can  be  studied  as  a  two-dimensional  problem,  one  dimension 
being  along  the  field,  i.e.,  along  z,  and  the  ether  transverse  to 


-  50  - 


the  field,  along  y.  Now  suppose  that  a  collision  has  just  taken 

place,  at  time  t^  and  at  the  co-ordinates  y  =  y^  and  z  =  ,  and 

that  we  have  determined  the  two  velocity  components,  v  .  and  v 

y,l  Z,1 

with  which  the  electron  emerges  from  the  collision.  The  computer 
now  generates  a  random  number  r,  and  (10)  then  detej  mines  the 
path  length  s  to  the  next  collision.  It  remains  to  determine  the 
electron  co-ordinates,  y  =  y^,  and  z  =  z^,  at  this  collision,  as 
well  as  the  electron  energy  E^.  In  order  to  do  this  it  is 
necessary  first  to  determine  the  time  t  of  the  second  collision. 
This  is  given  by  solution  of  the  integral  equation 


'  2 .  qj 

v  .  +  (v  +  t  ) 
y  ,1  z,i  m 


with  q  the  electronic  charge,  m  the  electronic  mass,  F  the  applied 
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field,  and  with  s  given  by  (10).  In  his  earlier  work,  Martinelli 
solved  this  equation  numerically,  by  the  technique  of  successive 
iteration.  [Noting  that  the  (indefinite)  integral  in  (ll)  is 
basically  of  the  form  [  dT  f  1  +  j2  =  i  T  ♦  L  ln  {T  ), 


we  see  that  the  integral  equation  (ll)  can  indeed  be  replaced  by 
an  algebraic,  transcendental  equation.  However,  this  must  obviously 
be  solved  by  some  method  of  successive  approximation  and  consequently 
there  may  not  be  any  substantive  saving  in  computer  time  over  Martinelli 's 
direct,  iterative  solution  of  (ll)].  Having  now  determined  t  ,  and 
setting  t«*  tf  -  t^,  the  electron  co-ordinates  and  energy  at  the  time 
tf  of  the  next  collision  are  given  respectively  by 


z  —  z  =  v  .t  +  — 
f  i  z,i  2m 


(12a) 


yf  '  yi  =  Vy,it 


(12b) 


E  «  j  m  [v  2  +  (v  +  t)2] 

1  ^  y,i  z,i  m 
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The  electron  loses  the  energy  through  single  phonon 
emission  at  each  collision,  and  then  starts  on  a  new  trajectory 
in  a  randomly  chosen  direction,  determined  by  a  computer-generated 
random  number  in  the  manner  discussed  below.  The  cycle  is  repeated 
until  the  electron  emerges  from  the  film  at  z  =  t.  An  energy 
(and  angular )  distribution  for  the  emerging  electrons  is  then 
obtained  by  repeating  this  entire  process  for  a  reasonably  large 
number  of  electrons,  generally  not  less  than  200  electrons,  and 
in  certain  cases  more. 

We  have  described  above  the  genera1  scheme  of  the  Monte 

Carlo  calculation  for  our  problem  (with  additional  details  supplied 

20 

below).  Martinelli's  earlier  experience  with  his  relatively 
simple  problem  made  it  clear  that  without  further  simplification 
a  comprehensive  computer  study  of  the  hot-electron  problem  would 
be  too  expensive.  The  simplification  that  we  decided  to  employ 
is  the  ingenious  self-scattering  technique  of  Rees  which  has 
proven  itself  eminently  useful  in  the  study  of  the  Gunn  effect 
and  related  phenomena  at  the  microscopic  level. 

If  there  is  a  constant  (energy  independent)  mean  free 
time  (mft)  T  between  successive  collisions  then,  by  direct  analogy 
with  (8b)  we  can  write  for  P(t),  the  probability  that  the  electron 
has  gone  for  time  t  without  scattering, 

P(t)  =  exp  (•  ^  .  (lUa) 

not  a  constant , 

In  this  case 

(lUc) 

which  is  the  direct  analogue  of  (8a). 

The  quantity  t()c)  is  normally  related  to  microscopic  processes 

through  standard  perturbation  theory  which  yields  the  relation: 


More  generally,  the  scattering  time  t  is 
but  depends  on  the  (crystal)  momentum  t  *t()j). 
P(t)  is  given  by  the  more  general  relation. 


P(t)  =  exp  - 


_dt 

T 
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wi  *  V&1  =  =  1  d*E’  sx(>S*)5')  •  (15) 

where  S^kjJc')  d)c'  is  the  transition  probability,  per  unit  time, 
for  the  electron  to  scatter  from  state  k  to  states  in  the  interval 
d}c'  about  Jc*  .  The  subscript  '1'  has  been  added  to  identify  the 
quantities  as  arising  from  a  given  scattering  process,  of  type  1. 

The  quantity  is  called  the  scattering  frequency  (for  processes 
of  type  l) . 

In  analogy  with  P  (t),  previously  defined,  we  define 
sc 

P  (t)  as  the  probability  density  that  the  electron  scatter  in 
s  c ,  1 

the  incremental  time  interval  dt  after  moving  for  the  time  t 
without  scattering.  In  analogy  with  (6a)  and  (6b)  we  have, 
respectively: 

P  (t)  /  t  \ 

p  (t)  =  _i -  a  WP(t)  =  W.  exp/-  W  (t')dt’J,  ( l6a) 

sc,l  11  1  V  J  o  1  / 

and,  for  constant  mft , 

p  (t)  «  —  exp  (-—-)=  W  exp  ( -W  t )  .  (l6b) 

sc,l  ^  ^  JL 
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If  there  are  different  scattering  processes  then, 
under  the  assumption  that  they  operate  independently  of  each 
ether,  the  total  scattering  frequency  W  is  the  sum  of  the  scatter¬ 
ing  frequencies  characterizing  the  separate  processes: 

W  =  ”  and  Eq.  19  generalizes  simply  to 


.  I,  W  , 

1=1  l 


r*  exp  - 


rt  n 

I  dt'  Z  W.(t') 
•3  i«l  1 


For  the  simple  constant  mft  problem  the  integration 
in  (19)  is  easily  performed  to  yield  the  direct  analogue  of  (10), 
namely 

t  =  — t  £n  r  •  (21) 

In  general,  of  course,  the  integration  in  (20)  cannot 
be  performed  analytically.  Here  the  ingenious  idea  of  Reas  is 
brought  into  play;  namely,  a  fictitious  self-scattering  process 
is  introduced  such  that  when  it  is  included  in  the  sum  I  VT  the 
integration  in  (20)  can  be  performed  to  yield  a  simple  result. 

The  choice  of  the  self-scattering  transition  probability 

is  straightforward: 


S  (k,k')  =  (r  -  I  W.)  <S(k  -  ic*  ) 
°  '  i=l  1 


where  r  is  an  arbitrary,  constant  frequency  which  need  only 

exceed  "  IT  and  the  delta  function  is  included  to  ensure 
' =1  Wi * 

self-scattering. 


The  corresponding  self-scattering  frequency  W  is 


given  by 


W  =  I  dk'  S  (k,k' )  =  T  -  E  W 


o  V  V 


But  (23)  can  be  re-written  in  the  form 


r  =•£«  w< 
1=0  l 
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If  now  the  integrand  in  (20)  is  extended  to  include 
the  self-scattering,  i.e,,  the  sum  "  w  is  replaced  by  "  , 

gives  the  simple  result 

t  =  -  j  in  r  (25) 

With  the  parameter  T  chosen  at  the  outset,  the  computer¬ 
generated  random  number  r  determines  the  interval  t  between 
collisions.  It  still  remains  to  determine  the  nature  of  the 
collision,  whether  self-scattering  or  physical.  This' requires 
another  random  number.  We  illustrate  this  for  the  simple  model 
of  a  single  physical  scattering  mechanism  which  we  assume 
corresponds  to  a  constant  mfp  X.  The  corresponding  scattering 
time  t  is  energy-dependent  and  is  given  by 


Ti  =  Ti(E)  =  7W 


(26) 


where  v(E)  denotes  the  electron  velocity  at  energy  E.  The 

choice  of  T  must  be  such  v.hat  T  >  — t=-t  =  W  (E)  for  most 

1 

energies  E  that  will  be  encountered.  A  new,  computer-generated 
random  number  r'  will  now  determine  whether  the  collision  is 


physical : 


0  <  r '  <' 


W.  (E) 


(27a) 


or  self-scattering:  W  (E )  <  <  ^ 

f 


(27b) 


If  the  collision  is  self-scattering  the  particle  proceeds 

in  its  Newtonian  trajectory  as  if  it'  had  not  made  a  collision 

(hence  the  terminology,  'self-scattering').  If  the  collision  is 

'physical'  it  is  allowed  provided  E(t)  >  e  ,  where  e  ,  is  the 

ph  ph 

assumed  optical  phonon  energy.  If  E(t)<  e  ,  •  the  collision' is 

ph 

treated  as  a  self-scattering  one.  For  the  usual  situation, 

E(t)>  ep^»  t*16  electron  emerges  from  the  collision  with  energy 
E(t)  -  e  k*  and  we  need  further  random  numbers  (as  discussed  below) 
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tc  determine  the  electron  trajectory  following  the  collision. 

Note  that  the  choice  of  a  suitable  T  is  not  critical. 

It  should  be  large  enough,  but  not  too  large.  If  it  is  too  large, 
an  inordinate  amount  of  computer  time  is  spent  generating  self¬ 
scattering  collisions,  which  contain  no  physical  information.  If 


cne  chooses  a  reasonable  r  which  then,  in  mid- calculation ,  proves 
tc  be  not  large  enough,  one  can  increase  its  value  to  a  new  r 

without  impairing  the  calculation. 

IT  remains  to  discuss  the  angular  distribution  of  the 

immediately  following  the  collision.  Talcing  the  k-vector  of  the 
electron  just  prior  to  the  collision  to  define  a  polar  axis,  the 
element  of  solid  angle  defining  the  post-collision  trajectory 
can  be  written  as 


d ft  =  sin9d9d$  =  dud$  »  ^  -  cos9  » 

where  9  is  the  polar  angle,  i.e-,  the  angle  between  Jc  and  k ,  and 

0  is  '“'he  azimuthal  angle- 

If  pf  u  J  is  the  normalized  probability  density  {in  fi)  fcr 
scattering  into  the  angle  u,  then  the  connection  with  a  computer¬ 
generated  random  number  r  in  the  interval  0  to  1  is: 


,  L  fl 

r  =  I  P(u)  du  ,  with  I  P(m)  dv  =  1-  (29) 

•'-1  '-1 

The  several  cases  we  have  studied  are  as  follows: 

( i )  Isotropic  scattering 

Here  all  u  are  weighted  equally  so  that  (29)  gives 
p  (u)  =  j  ,  independent  of  V ,  r  =  77  (  P  +  !)»  (30a) 


cr 

y  =  2r  -  1. 


(30b) 


The  subscript  I  in  (30a)  is  used  to  denote  isotropic  scattering. 
Likewise  all  <t>  are  weighted  equally  so  that 
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PU) 


_1 

2tt 


P(*)d*  =  ^ 


(31a) 


or 


<P  =  2nr  .  (31b) 

(ii)  Dominant  Forward  Scattering 

Here  we  have  arbitrarily  taken  the  simplest  forward 
scattering  law,  namely  Pp(cos0  )  «c  1  +  cos  0  ,  to  represent  sche¬ 
matically  the  effects  of  dominant  forward  scattering: 

pp(u)  =  j  (l  +  y)  ,  r  =  j  y2  +  y  +  )  ,  (32a) 

or,  solving  for  y  : 


y  =  2  /r  -  1  .  (32b) 


The  distribution  of  the  azimuthal  angle  <j>  is  again  given 
by  (31a)  and  (31b). 
iii)  Dominant  Backward  Scattering 

As  with  forward  scattering,  we  have  taken  the  simplest 
backward  scattering  law  to  represent  schematically  the  effects  of 
dominant  backward  scattering,  namely  P_(cos0)  «  1  -  cos  0,  or 


P b(m)  ■  \  (l  -  w) 


1  , 12 
r  =  2  (2  y 


(33a) 


or,  solving  for  y  , 


y  =  1  -  2^7  .  (33b) 

The  distribution  of  the  azimuthal  angle  $  is  again  given 
by  (31a)  and  (31b). 
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1V)  Polar  Optical  Phonon  Scattering 

Here  the  starting  point  is  the  scattering  relationship 


WP0( '  '  Ik'-k 


const 

r 


const _ _ 

k'^  +  k'  -  2kk'u 


(34) 


Through  some  straightforward  algebra  we  obtain  the 
desired  relation  of  u  with  the  random  number  r  : 


In 

r  (k  +  k')?  1 

1  fk-k' )2  +  2kk ' ( 1-  y)  . 

II 

H- ' 

D 

f  (k  +  k')  2V 

C\J 

_  / 

Conservation  of  energy,  for  phonon  emission,  is: 


E  = 


2  2 

h  k 

2m 


2  2 

ft  k' 

2m 


+ 


(36) 


Further  a  gebra  enables  us  to  express  conveniently  y 
in  terms  of  r: 

v2 


y  ■  1  -  “ 


SrT  +  /r-l)^  exp(  fr. )  -  (  v/iT  - 


2  /n  'Ai  -1 


,  (37a) 


with 


r 


n 


r  In 


/  y/ri~~  ■f'/r-l 
(  '/p"  -/r-1 


) 


2 


(37b) 


C.  Computer  Results  for  Constant  MFP  Problems 

(i )  Three-Dimensional  Computations,  Isotropic  Scattering 
The  above  discussion  presents  in  considerable  detail  the 
basic  ingredients  that  have  gone  into  our  Monte  Carlo  studies  of 
the  hot-electron  problem.  Beyond  this  point  we  are  dealing  with 
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straight  programming  which  need  not  be  discussed  here.  We  there¬ 
fore  proceed  directly  to  the  results.  In  Fig.  15  is  plotted  the 
average  energy  E  of  the  electron  emerging  from  the  film  as  a 
function  of  film  thickness  t,  for  the  situation  of  constant  mfp 
and  isotropic  scattering.  Six  different  cases  are  plotted, 
corresponding  to  two  different  starting  energies:  E  =  O.Olev  and 

O  £- 

Eq  =  l6  ev,  three  different  field  strengths:  F  =  5  x  l(r  V/cm, 

F  =  1  x  10  V/cm  and  F  =  2  x  10  V/cm,  and  three  different  mfp's: 

0  0  0 

^  =  25A,  A  =  50A  and  A  =  100A.  However,  for  all  cases  FA  is 

the  same ,  namely  FA  =  0,5  V;  likewise  e  ,  is  common  to  all  cases 

Ph 

namely  e  =  0.1  ev.  It  is  seen  that  ,  at  sufficiently  large 

>  O 

thicknesses  t  ^  U500A^all  of  the  cases  have  converged  to  the 
common  value,  E&ve  sg=  2.1  ev.  (subscript  'ss'  denotes  steady 
state).  This  one  plot  already  leads  to  the  conclusion  that 

n  pp 

E  a(FA)  .  From  theoretical  work  of  Wolff  ,  as  well  as 

clVG  9  S  S  * 

further  computer  studies  (see  Figs.  23,  20,  2.0),  the  more  precise 
relation  follows : 


> 


E 

ave ,ss 


=  k 


k  =  0.87 


(38) 


This  relation  had  already  been  found  in  the  earlier  computer 

20 

studies  of  Martinelli  ,  who  determined  k  =  0.88.  We  have 

shown  (38)  to  hold  (with  a  variance  of  k  less  than  10$)  over 

the  range  of  energies  0.5  ^  E  v  18  ev.  (An  actual 

2  & v e  ,ss 

plot  of  Eave  ss  vs.  (F  A)  /£^h  is  given  in  Fig.  23,  the  middle 
straight  line). 

It  has  also  been  found  through  the  computer  studies 
that,  for  the  constant  mfp  isotropic-scattering  problem,  the 
electrons  warm  up  or  cool  down  to  the  steady-state  energy  in  an 
exponential  manner: 


E  -  E  *  (E  -  E  ) 

ave  ave  ,  s  s  0  ave , s  s 


(39) 


AVERA6E  ENERGY  OF  EMERGENCE 
ISOTROPIC  SCATTERING  . 


NUMBER  OF  PARTICLES  •  500 
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Vith  Eo  the  initial  energy  of  the  electrons,  at  the  entrance  plane 
t  =  0. 

In  our  computer  studies  to-date  we  have  taken  E  as  mono- 
energetic;  if  the  electrons  enter  with  an  energy  distribution,  then 
presumably  (39)  still  holds  (at  least,  approximately)  with  E  re¬ 
placed  by  The  development  distance  obeys  an  approximate 

scaling  law  similar  to  (38),  namely 


D  =  k' 


with  1  <  k'  <  2.  (1*0) 


20 

Martinelli  ,  in  his  more  restricted  study,  found  k'  =  1.27.  The 
development  distance  for  the  six  cases  considered  in  Fig.  19  is 
established  by  the  plots  in  Fig.  20. 

(ii)  Three  Dimensional  Computations.  Anisotropic 
Scattering 


In  addition  to  the  isotropic  scattering  case,  anisotropic 

scattering  cases  have  also  been  studied  for  the  constant  mfp 

problem.  In  Fig.  21  is  plotted  the  average  energy  E  of  the 

ave 

emerging  electrons  vs.  film  thickness  t  for  two  different  examples 
of  forward  scattering:  moderate  forward  scattering  characterized 
by  a  (1  +  cos  9)  distribution  and  with  y  related  to  r  through 
(32a)  and  (32b),  polar-optical  forward  scattering  as  given  by 
(3*0  thru  (37b).  For  each  of  the  two  examples  we  have  taken 


two  widely  separated  starting  energies,  a  low  energy  and  a 

high  energy.  We  see  that  there  is  indeed  a  steady-state  average 

energy  for  each  of  the  two  examples,  and  that  this  energy  is 

higher  for  the  polar-optical  case,  Eo_  -  12.2  ev,  than  for 

the  (1  +  cos  9)  case,  E^^  =  3.0  ev,  both  being  calculated 

for  the  same  parameter  values:  \  =  50A,  e  =  0.1  ev  and 
6  ph 

F  -  1  x  10  V/cm.  This  result  is  expected  since  polar-optical 

scattering  is  more  strongly  peaked  in  the  forward  direction.  The 

vertical  line  bracketing  each  average-energy  point  is  not  an  error 

bar;  rather,  it  gives  directly  the  half-width  of  the 


energy 
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distribution,  [in  a  precise  formulation  of  the  polar-op  ical 

scattering  problem  the  mean  free  path  is  not  a  constant,  but  is 

energy-dependent,  increasing  with  increasing  energy.  We  have  an 

approximate  formulation  of  energy-dependent  mean-free-path  problems, 

discussed  below,  which  enables  us  to  study  them  first  as  constant 

mean-free-path  problems.]  A  corresponding  plot  for  moderate 

backward  scattering  characerized  by  a  (l  -  cos  9)  distribution 

and  with  related  to  r  through  (33a)  and  (33b)  is  presented  in 

Fie  22.  The  >  e  and  F  are  the  same  as  in  Fig.  21.  The  value 

&  ph 

of  E  is  now  1-6  ev,  lower  than  the  two  values  for  the 

ave ,ss 

forward  scattering  examples  of  Fig.  21,  as  expected. 

We  have  already  noted  the  proper  scaling  law  for 

E  ,  for  the  isotropic  scattering  case,  in  Eq.  (38).  The 

ave  ,ss 

same  scaling  law  holds  for  the  moderate  anisotropies :  forward 
scattering  characterized  by  a  (l  +  cos6  )  distribution,  and 


backward  scattering  characterized  by  a  (1  -  cos  0)  distribution, 

as  seen  by  the  straight-line  plots  for  the  three  kinds  of 

scatterings  in  Fig.  2£.  The  corresponding  proportionality 

constants  kT ,  kn  (I  for  isotropic,  F  for  forward  and  B  for 

backward)  are  KT  =  0  87,  k  =  1 .27  and  k  -  0. 66 .  Note  that 
I  F  D 

k  <  kT  k  ,  as  physically  expected.  An  approximate  scaling 
B  I  F 

law  exists  for  pclar  optical  scattering  (however,  with  constant 
mfp),  and  is  mentioned  below. 

The  fact  that  anisotropic  scattering,  no  matter  how 
strongly  forward  [short  of  pure  (100JS)  forward  scattering]  by 
itself  (that  is,  without  an  energy- depen dent  mfp)  leads  to  a 
stable  distribution  has  been  verified  by  calculations  witv  a 
one-dimens  onal  model,  as  described  below.  This  situation  can 
be  understood  theoretically  by  noting  that  even  though  the 
probability  for  forward  scattering  along  a  trajectory  is  large, 
there  is  a  finite  probability  that  the  electron  be  backward 
scattered.  Once  a  high  energy  electron,  the  trajectory  of  which 
is  not  appreciably  affected  by  the  field,  is  back-scattered  in 
the  direction  against  the  field,  it  will  continue  to  travel  in 
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that  direction  (as  the  electron  has  high  probability  of  being 
scattered  forward  along  the  trajectory)  until  the  small- 
probability  event  that  it  be  scattered  in  the  opposite  direction 
from  the  direction  of  the  trajectory  occurs,  or  until  it  loses 
enough  energy  such  that  the  field  can  turn  it  around  along  its 
own  direction.  The  existence  of  this  cooling  mechanism 
stabilizes  the  electron  energy  distribution  and  a  steady 
state  thus  exists. 

(iii)  One-Dimensional  Simulations 

Even  with  the  saving  in  computer  time  gained  through  use 
of  the  self-scattering  technique, the  three-dimensional  calculations 
required  with  anisotropic  scattering  are  still  quite  costly. 
Therefore  further  simplification  is  highly  desirable.  With  this 
aim  in  mind  we  decided  to  explore  the  usefulness  of  one-dimensional 
models  in  simulating  the  full  three-dimensional  problems  .  The 
simplest  possible  one-dimensional  model  is  the  random  walk  in  one 
dimension.  The  particle  is  constrained  to  move  on  a  line  and  after 
moving  a  fixed  distance  A (the  same  A  which  characterizes  the  con¬ 
stant  mfp  problem  in  three  dimensions)  along  its  trajectory  on 
the  line  it  makes  a  collision.  Note  however  that  this  is  not  the 
classical  random-walk  problem,  which  deals  with  force-free 
situations.  Here  the  electron  is  acted  on  by  the  eleetric  field 
(directed  along  the  line)  throughout  its  motion.  In  other 
words,  our  problem  is  a  random  walk  in  a  force  field.  The  first 
problem  of  this  kind  which  we  studied  is  the  analog  of  the 
isotropic  scattering  problem.  Here,  and  in  the  following,  we 
shall  always  be  referring  to  constant  mfp  situations,  both  in 
three  and  in  one  dimensions.  The  obvious  one-dimensional  analog 
of  isotropic  scattering  is  equal  probabilities  for  scattering 
forward  and  backward  n  the  line  respectively.  In  Fig.  24  we 
compare  the  spatial  evolution  of  the  average  energy  E 

the  two  problems  with  the  same  paramete-s:  E  =  0.01  ev, 
o  /’  ° 

A  =  50A  ,  e  ,  =  0.1  ev,  and  F  =  1  x  10  V/cm. 

ph 
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-i.n  Fig.  25  the  spatial  evolutions  of  E  for  the  two  problems 

3.  V  G 

are  compared  for  the  same  parameters  except  for  E  which  is 

0 

now  a  higher  energy.  Figures  24  and  25  indicate  how  good  a 
simulation  the  one-dimensional  random  walk  furnishes  for  the 
three-dimensional  problem.  The  one— dimensional  curves  in 
Figs.  24  and  25  are  re-plotted  in  Fig.  26  to  show  their 
convergence  to  a  common  value  of  E  .  The  development 

distance  for  five  cases,  in  one  dimension,  is  established  by 
the  plots  in  Fig.  2 7  (the  analog  to  Fig.  20). 

The  problem  of  finding  the  most  suitable  one¬ 
dimensional  random  walk  (in  a  force  field)  to  simulate  a 
given  anisotropic  three-dimensional  problem  is  one  step 
more  subtle.  The  added  complexity  comes  from  the  need  to 
establish  the  proper  value  for  f,  the  probability  for 
'forward'  scattering  on  the  line  (1-f  is  then  the  probability 
for  'backward'  scattering  on  the  line).  After  trying  unsuc¬ 
cessfully  to  establish  the  proper  value  of  f  by  guessing 
'appropriate'  averages  for  the  three-dimensional  anisotropy, 
we  switched  to  a  very  simple,  straightforward  trial-and- 
error  method.  We  have  already  noted,  in  Fig.  23,  that  for 
each  of  the  three-dimensional  problems:  isotropic,  (l  +  cos  0) 

and  (l  -  cos  0)  scattering, E  «(F  \)2 jz  ,  .  It  would 

f  ave j  s  s  pn 

seem  reasonable  that ,  whatever  the  value  of  f  between  0  and  1 

the  same  scaling  law  would  hold;  but  with  a  different 

constant-of-proportionality  k^  corresponding  to  each  f.  That 

this  is  indeed  true  is  seen  in  Fig.  28.  Knowing  the  k  of  the 

scaling  law  (assumed  given  by  Eq.  (38)  for  any  given  anisotropic, 

three-dimensional  problem  it  is  only  necessary  to  find  that  f 

for  which  kf  is  the  same  k.  Using  this  criterion,  we  find  from 

Fig.  28  that  the  appropriate  f  for  the  one-dimensional  simulation 

of  the  (l  +  cos  0 )-anisotropy  is  f  =  0.6,  and  that  the  appropriate 

f  for  the  simulation  of  the  (l  -  cos  0)-anisotropy  is  f  =  0.^. 

V/e  would  also  expect  that  the  spatial  development  of  E  to  its 

ave 

steady-state  value  Eavg  ss  is  given  by  Eq.  (39),  with  the  same 
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scaling  law  (40)  for  D,  only  a  new  constant  k^  corresponding  to 

the  new  value  of  f  That  this  is  indeed  true  is  seen  in  Fig.  29 

2 

which  gives  a  plot  of  the  development  distance  D  vs.  F  A  /e 
for  the  three  one -dimensional  cases:  f  =  0.6  (moderately  forward 
scattering)  and  k'  =  1-40,  f  =  0,5  (isotropic  scattering)  and 
k'  =  0.94  [  compared  to  1  -  k'  2  for  the  three-dimensional 
isotropic  case],  and  f  =  0.4  (moderately  backward  scattering) 
and  k'  =0.68. 


A  direct  comparison  of  three-dimensional  (l  +  cos  0)- 
scattering  and  f  =  0,6  one-dimensional  scattering  is  given  in 


Fig.  305  which  presents  a  plot  of  E&ve  vs 


film  thickness  t  for 


the  two  problems  under  'identical'  conditions:  E  =  0.01  ev, 
o  ^  o 

A  =  100A,  e  ,  =  0  1  ev,  and  F  -  1  x  10°  V/cm,  The  agreement 

ph 

is  seen  to  be  excellent-  Note  also  that  the  relative  widths 

of  the  two  distributions  is  almost  identical:  (  a/E  )  =0.53 

ave  ss 

and  0  51  for  the  three-and  one-dimensional  cases  respectively. 


Here 


C  "  /(  E2)  -  E 


and  ss ,  as  usual,  denotes  the 


ave 


ave 


steady  state  Figure  31  describes  the  approach  to  the  steady 

state  as  a  funcetion  of  film  thickness  for  five  cases  for  the 

one-dimensional,  f  =  0.6  problem. 

A  direct  comparison  of  three-dimensional  (l  -  cos  0)- 

scattenng  and  f  =  0  4  one-dimensional  scattering  is  given  in 

Fig.  32  for  the  two  problems  under  'identical'  conditions  for 

o 

two  choices  of  parameter  values:  E  =0,01  ev,  A  =  100A, 

/•  o  o 

e  ,  =  0.075  ev,  and  F  -  1  x  10  V/cm;  E  =  0.01  ev ,  A  =  50A, 
ph  ° 

=  0  1  ev.  and  F  =  1  x  10°  V/cm.  Again,  the  agreement  is 
ph 

seen  to  be  excellent.  Figure  33  describes  the  approach  to 
the  steady  state  with  film  thickness  for  five  cases  for  the 
one-dimensional,  f  =  0  U  problem. 

Before  leaving  the  constant  mfp  problem  we  show  two 
sets  of  results  for  the  polar  optical  phonon  scattering  problem 
treated  as  a  constant  mfp  problem  (  A  =  const,  but  the  anisotropy 
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of  the  scattering  correctly  described  by  (3^)  thru  (37b). 

I:,  Fig.  34-,  we  have  clotted  E  vs.  film  thickness  t  for 

"  ave  o 

three  cases:  A  -  12, 5A  giving  E  =  0.32  ev,  A  -  25A 

ave  )|o 

giving  E  =  1-89  ev,  and  A  =  kOA  giving  E  =  6.66  ev 

ave  ,ss 

(for  ail  three  cases,  Eq  =  0.01  ev,  e  ^  =  0.1  ev  and 
F  =  1  x  10^  V/cm),  These  calculations  yeild  the  scaling  law: 

E  =  K  ( F5  )8/3/e  ,5'3  with  k  -  3-5.  Figure  35 

ave.ss  po  Fh  po 

describes  the  approach  to  the  steady  state  with  film  thick¬ 
ness  for  five  cases  of  poiar -optical  scattering  with 
(arbitrarily  assumed)  constant  mfp-  Note  that,  unlike  the 
previous  problems  plotted  in  this  manner,  these  curves  are 
not  straight  lines  on  the  semi -log  plot. 

D.  Computer  and  Graphical  Results  for  Energy -Dependent 
MFP  Problems 

The  final  major  topic  to  be  considered  in  this  interim 
report  on  our  hot -electron  calculations  is  the  energy-dependent 
mean  free  path:  A  =  *. (E  Tics  is  a  particularly  important 
problem  because  the  origin  cf  instability  of  the  distribution 
function  (runaway)  lies  precisely  in  a  A(E)  which  increases 
with  E  sufficiently  rapidly  In  the  framework  of  the  self¬ 
scattering  technique  there  is  very  little  change  in  the 
calculation  as  already  described  for  constant  mfp.  Once  again 
the  time  interval  between  scatterings  is  given  by  (25).  We  now 
need  only  know  the  particle  energy  E  at  the  time  of  the  collision 
to  assert  that  the  relative  probabilities  for  real  scattering 
and  self-scattering  respectively  are  v(E)/Ta  (E;  and 
1-  (v(E)/I '  i E ) j  A  random  number  is  used  in  the  usual 
manner  to  determine  which  kind  of  scattering  has  occurred. 

Up  to  the  present  we  have  studied  two  simp’e  'proto¬ 
typical'  cases  of  energy-dependent  mfp:  1)  A  =  +  A  , 

and  2)  A  =  A '  +  *'  E/e  where  A  ,  A'  are  constant  mfp's  and 
,  o  o  ph  o  o 

A  ,  A  are  suitable  constants  with  the  dimension  of  length  (we 
11  o 

actually  used  A^  =  A^  =  1A).  Three  kinds  of  situations  have 
been  encountered:  l)  stability,  ii)  quasi-stability  and  iii)  run- 
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away.  Further  we  have  found  that  these  three  situations  can  be 
predicted  through  graphical  solution  of  the  simultaneous  equations 


*  =  '(E)  and  E  =  E  (  A F,  e  ) ,  (4i) 

ave,ss  ss  ph 

where  m  the  latter  equation  X  is  taken  as  a  constant . 

For  the  cases  studied  tc-date  (excepting  polar  optical 
scattering)  we  have,  more  specifically,  to  deal  with  the  pair: 

( FX 

*  =  /(E)  and  E  =  k  — .  (42) 

ave,ss  e 


If  there  is  only  one  solution  to  the  pair  of  simultaneous  equations 

then  there  is  a  stable  steady  state.  Such  is  the  situation  in 

Fig.  36,  where  we  have  taken  X  =  \  +  2/  E/  e~T  (angstroms).  The 

o  ph 

^'responding  E  vs  t  plots,  as  determined  directly  from  the 

self-scattering  calculations,  are  exhibited  in  Fig.  37-  It  is 

seen  that  the  calculated  E  values  are  very  close  to  those 

ave,ss 

predicted  by  the  graphical  intersection  points  in  Fig.  56> 

namely  the  heavy  dots  in  Fig.  36.  (The  two  dots  in  Fig.  3fc  for 

o 

the  Xq  =  25A  case  correspond  to  the  different  starting  energies  in 

Fig.  37:  E  =  4.0  ev.  and  E  =  0.01  ev.  The  two  dots  should 
o  o 

coincide;  their  difference  is  small  enough  to  lack  significance. 

Figure  36  indicates  a  situation  where  there  are  two 

intersections  or  none,  namely  X  =  X  +  X.  E/e  .  and  there  are 

o  1  ph 

two  for  1  and  none  for  A  >1.  In  the  two-intersection  cases, 

the  two  points  of  intersection  are  wider  apart  the  smaller  is 
For  near  1  the  intersection  points  are  close  together  and 
the  corresponding  Eq  ”s .  t  plots,  exhibited  in  Fig.  39  are  quite 
erra+ic  -  as  if  the  distribution  has  difficulty  making  up  its 
mind  on  what  to  do  -  but  appears  headed  for  runaway.  This  unusual 
situation  (two  close  intersection  points)  certainly  needs  further 
study.  On  the  other  hand  with  widely  separated  intersection 
points,  X  =  0.5,  we  see  from  the  Eavg  vs.  t  plots  in  Fig.  40 
that  there  is  at  least  quasi-stability.  Finally,  with  X^  =  2  in 
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Fig,  3S»  there  are  no  intersection  points  and  accordingly  the 
distribution  simply  runs  away,  as  is  clearly  visible  in  the 
E  vs-  t  Plot  in  Fig.  41. 

The  energy-dependent  mean-free-path  studies  reported 
here  are  fragmentary  in  nature.  Since  this  case  is  technologi¬ 
cally  an  important  one,  more  systematic  studies  will  be  made. 

It  remains  to  be  seen  whether  the  one-dimensional  simulations, 
either  random  walk  or  non-random  walk,  will  prove  as  useful 
for  the  study  of  this  class  of  problems 
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